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PREFACE 


In  January  1985,  Dr.  James  A.  Ionson,  then  Director  of  the  Science  and 
Technology  Directorate  of  the  Strategic  Defense  Initiative  Organization  (SDIO)  asked 
members  of  the  IDA  research  staff  to  investigate  the  feasibility  of  developing  a  gamma-ray 
laser.  The  staff  determined  what  work  had  been  done,  who  was  currently  working  in  the 
field,  and  what  work  should  be  encouraged  or  supported.  Then  a  workshop  was  convened 
for  research  workers  directly  involved  both  in  gamma-ray  laser  work  and  in  ancillary  fields 
such  as  nuclear  structure,  radiation  propagation  in  crystals,  Mossbauer  effect,  and  optical 
lasers.  Next,  an  in-house  study  was  undertaken  to  clarify  critical  issues  concerning  the 
various  pumping  schemes  proposed  at  the  workshop  as  well  as  systems  questions  about 
the  gamma-ray  laser  as  a  working  device. 

The  proceedings  of  the  workshop  were  published  in  the  form  of  a  report  to  the 
Innovative  Science  and  Technology  Office  (1ST)  of  the  SDIO.  The  work  completed  in 
1985  is  presented  in  IDA  Paper  P-2021,  "IDA  Gamma-Ray  Laser  Annual  Summary  Report 
(1985):  Investigation  of  the  Feasibility  of  Developing  a  Laser  Using  Nuclear  Transitions" 
(Ref.  1). 

In  1986,  the  in-house  work  focused  on  extending  the  data  base,  the  nature  of 
superradiance  in  the  gamma-ray  laser  context,  and  a  detailed  investigation  of  the 
upconversion  pumping  scheme.  A  discussion  of  nuclear  systematics,  investigations  of 
electron-nuclear-driven  pumping,  and  lifetime  measurements  rounded  out  that  study.  The 
results  of  the  FY  1986  effort  are  presented  in  IDA  paper,  P-2004  (Ref.  2). 

The  development  of  a  gamma-ray  laser  is  viewed  as  a  high-risk/high-payoff 
undertaking.  IDA's  involvement  focuses  on  minimizing  the  risk  and  on  striving  to  redirect 
the  effort  when  proposed  schemes  are  shown  not  to  be  feasible. 

In  1987,  the  IDA  research  staff  looked  at  the  state  of  the  art  and  assessed  the 
situation  in  gamma-ray  laser  work  (described  in  Chapter  I).  Focus  was  on  two  areas  of 
research  interest  critical  to  three  concepts  for  developing  a  gamma-ray  laser.  In  Chapter  II 
the  heating  effects  associated  with  upconversion  techniques  are  discussed.  Although  the 
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work  concentrated  on  the  direct  upconversion  of  nuclei  with  electromagnetic  radiation,  the 
results  of  the  studies  also  apply  to  upconversion  through  the  mechanism  of  electron-nuclear 
coupling.  In  the  second  study  (Chapter  HI),  the  sources  of  inhomogeneous  broadening 
which  destroy  the  Mossbauer  effect  are  investigated  and  the  techniques  available  for 
restoring  the  resonance  (destroying  the  effect  of  inhomogeneous  broadening)  by  external  or 
internal  fields  are  considered.  An  old  concept  for  doing  this  is  reviewed  and  two  new 
approaches  are  introduced. 

This  report  does  not  have  an  overall  introduction.  Each  chapter  is  an  independent 
study  containing  its  own  introduction. 
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SUMMARY 


This  report  describes  the  1987  research  effort  by  members  of  the  EDA  staff  in  the 
field  of  gamma-ray  lasers.  The  work  is  part  of  a  continuing  task  in  support  of  the 
Innovative  Science  and  Technology  Office  (1ST)  of  the  Strategic  Defense  Initiative 
Organization  (SDIO).  The  development  of  a  gamma-ray  laser  is  a  high-risk  science  and 
technology  undertaking.  IDA  involvement  has  focused  in  large  measure  on  minimizing  the 
risk  and  attempting  to  redirect  the  program  as  quickly  as  possible  when  proposed  schemes 
prove  infeasible.  The  report  is  presented  in  three  independent  chapters. 

Chapter  I  is  a  tutorial  overview  of  the  subject  of  gamma-ray  lasers.  It  presents  a 
review  of  the  problems  specific  to  nuclear  transitions  that  make  the  task  of  building  a 
nuclear  laser  more  difficult  than  building  conventional  lasers.  A  unique  two-step  approach 
is  taken,  showing  first  that  a  gas  laser  is  possible,  in  principle,  then  considering  a 
crystalline  laser.  Then,  assuming  that  a  gamma-ray  laser  is  feasible,  major  problem  areas 
requiring  study  are  pointed  out  These  include  establishment  of  laser-specific  nuclear  data 
bases,  production  of  samples  of  the  right  nuclear  species,  fabrication  of  the  required 
crystals,  and  detailed  studies  of  radiation  transport  and  coherent  emission.  The  three  major 
lasing  schemes  that  have  been  proposed  are  discussed,  along  with  the  specific  problems 
that  will  have  to  be  overcome  in  each  case. 

Chapter  II  takes  a  critical  look  at  the  heating  effects  associated  with  upconversion 
techniques.  Sample  calculations  for  x-ray  upconversion  in  very  thin  wires  of  57Fe  lead  to 
phase  changes  from  heating  which  destroy  the  lattice  and  therefore  the  Mossbauer  effect. 
Preliminary  calculations  indicate  that  conductive  cooling  of  the  wire  is  not  a  viable  option. 
It  is  therefore  suggested  that  other  pumping  schemes  receive  greater  attention. 

Chapter  HI  reexamines  the  long-lifetime  concept  of  developing  a  gamma-ray  laser. 
The  major  difficulty  with  this  scheme  is  inhomogeneous  broadening.  The  problem  is 
approached  from  two  extremes.  In  one,  the  broadening  is,  in  principle,  removed  through 
the  use  of  RF  pulses  which  undo  the  effects  of  the  dipole-dipole  interactions  in  the  crystal. 
In  the  other  approach,  the  level  is  homogeneously  broadened  through  level  mixing 
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spectroscopy  techniques  or  through  relaxation  effects.  The  enhanced  broadening  increases 
the  overlap  of  transitions  and  may  produce  an  increase  in  the  nuclear  resonance.  Further 
work  has  to  be  done  to  determine  the  effect  on  the  nuclear  cooperative  emission  or 
supeiradiance. 
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I.  ASPECTS  OF  THE  CURRENT  STATUS  OF  THE  GRASER 
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A.  INTRODUCTION 

At  the  core  of  every  atom  is  a  small  entity,  the  nucleus.  The  radius  of  the  nucleus  is 
of  the  order  of  a  few  femtometers,  i.e.,  it  is  about  five  orders  of  magnitude  smaller  than  the 
atom.  It  is  a  system  of  many  particles  and  is  amenable  to  a  quantum-mechanical 
description.  The  particles  give  rise  to  charge  distributions  and  currents  which  are  the 
sources  of  the  electromagnetic  radiation  that  we  are  interested  in.  This  radiation,  both  in  its 
multipole  character  and  in  its  interactions  with  matter,  is  quite  different  from  the  longer 
wavelength  radiation  emitted  during  changes  in  the  electronic  structure  of  the  atom.  As  in 
the  case  of  the  atom,  the  nucleus  is  characterized  by  a  complex  eigenvalue  spectrum.  Our 
goal  is  to  study  the  feasibility  of  a  laser  based  on  an  electromagnetic  transition  between  two 
of  the  levels  in  a  particularly  suitable  nucleus. 

Figure  1  shows  a  modification  of  an  historical  plot  by  D.  J.  Nagle.  (Ref.  3)  The 
graph  shows  the  peak  energy,  as  a  function  of  time,  of  new  sources  of  coherent  radiation. 
The  two  highest  points  (represented  by  squares)  and  a  shaded  horizontal  band  have  been 
added.  The  points  represent  two  reported  x-ray  lasers  and  the  band  shows  the  approximate 
region  of  interest  for  the  prospective  graser,  roughly  from  less  than  100  eV  to  more  than 
100  keV.  A  major  reason  for  the  interest  in  nuclear  lasers  is  the  availability  of  a  region  of 
higher  photon  energy  and  the  concomitant  possibilities  of  higher  densities  of  energy  storage 
and  more  intense  and  more  penetrating  beams.  To  see  the  contrast  between  nuclear  and 
atomic  structure,  it  is  worth  observing  the  level  structure  of  two  typical  nuclei,  one  light 
and  one  heavy.  Figure  2  is  an  energy  level  diagram  of  llC.  The  contrast  to  an  atomic  case 
is  rnmediately  obvious  as  excitation  energies  are  of  the  order  of  MeV,  not  eV.  The  lower 
leveis  decay  only  by  photon  emission,  but  the  higher  states  have  other  modes.  For 
example,  a  state  above  7.55  MeV  can  decay  by  the  emission  of  an  alpha  particle,  leaving  a 
7Be  nucleus.  Conversely,  these  states  can  be  formed  by  bombarding  a  7Be  target  with 
sufficiently  energetic  alpha  particles.  Figure  3  shows  the  energy  level  diagram  of  a  very 
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E(MeV) 


Figure  3.  The  level  spectrum  of  235U  is  shown  for  an  excitation  range  of  2 
MeV.  The  location  of  the  lowest  state  at  73  eV  is  exaggerated. 
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heavy  nucleus,  235U.  In  the  2  MeV  of  excitation  for  which  only  one  level  appeared  in  the 
carbon  nucleus  (Fig.  2),  there  are  now  about  100  levels  (Fig.  3);  and  the  first  excited  state 
has  an  energy  of  only  73  eV,  the  lowest  known  first  excited  state  of  a  nucleus.  This  level 
plays  a  role  in  one  of  the  schemes.  At  the  current  level  of  nuclear  theory,  these  spectra 
cannot  be  predicted  with  any  precision. 

Again,  our  goal  is  to  study  the  gamma-ray  laser— a  conceptual  device  based  either 
on  stimulation  or  superradiance.  It  is  based  on  a  transition  between  two  nuclear  levels.  As 
a  first  guess,  the  nuclei  will  be  embedded  in  a  crystal  or  gas  having  an  acicular  shape. 

B  .  MULTIPOLE  RADIATION 

The  sketch  in  Fig.  4  represents  radiation  of  frequency  to  that  is  emitted  when  a 
transition  occurs  between  the  two  levels. 


Eq  =hto0  =E2  -Et 

Figure  4.  The  two-level  system. 

The  power  that  is  radiated  can  be  expressed  as  a  series  of  two  terms: 


The  first  terms,  the  electric  multipoles,  are  due  to  the  charge  distributions  in  the  nucleus: 


The  second  set  of  terms,  the  magnetic  multipoles,  are  due  to  currents  and  are  given  by: 

Mt«  -  <«•♦)  nrf1  ix  -  ttt  m:  £ H  *1-  (W  ?  <v*2w<*. 


where 

=  -irk  x  Vk  .  (3) 


Each  term  gives  the  power  radiated  by  that  particular  electric  or  magnetic  multipole,  and  the 
corresponding  decay  constants  are  given  by 


V<E)  - 


pto<E) 

Ho) 


and 


P/m(M) 


(4) 


The  corresponding  lifetimes  are  given  by 


x 


l’m  X. 


1_ 

l,m 


(5) 


In  general,  the  number  of  terms  in  the  expansion  are  limited  to  just  a  few  by  selection  rules 
on  angular  momentum  and  parity.  These  are: 


Angular  Momentum:  |  J2-  J 1 1  £  L  ^  |  J2+ J 1 1 

Parity  Electric:  7C2Jti  =  (-1) L 
Magnetic:  7C2  7ti  =  (-1)  L+1 

Table  1  lists  the  L- values  and  indicates  whether  there  is  a  parity  change  for  the  first 
five  multipoles.  As  an  example,  if  there  is  a  transition  between  an  upper  state,  hn  =  1/2' 
and  Jin  =  7/2+,  the  possible  L- values  between  the  sum  and  difference  of  the  two  J’s  are 
just  3,4  ;  and  there  is  a  parity  change  (Yes).  Therefore,  according  to  the  table,  the  only 
possible  terms  in  the  multipole  expansion  are  E3  and  M4.  A  more  thorough  treatment  of 
the  multipole  expansion  can  be  found  in  standard  texts.  See,  for  example,  Ref.  4. 
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Table  1.  Parity  Changes  for  Multipole  Transitions 


El 

Ml 

E2 

M2 

E3 

M3 

E4 

M4 

E5 

M5 

D 

1 

2 

3 

4 

5 

Ait 

Yes 

No 

No 

Yes 

Yes 

No 

No 

Yes 

Yes 

No 

High  multipole  transitions  occur  frequently  in  the  nuclear  case  but  not  in  the  atomic 
case.  This  situation  arises,  in  part,  because  of  the  differences  in  atomic  and  nuclear 
structure.  In  the  atom,  electrons  find  themselves  in  the  long-range  coulomb  field  of  the 
positive  nucleus.  The  resulting  eigenvalue  spectrum  has  single  particle  levels  which  are 
filled  according  to  the  Pauli  exclusion  principle.  This  procedure  gives  rise  to  particularly 
stable  configurations  of  2,  8,  18,  etc.,  electrons.  In  the  nucleus,  there  are  two  types  of 
particles  which  are  different  charge  states  of  the  same  nucleon.  The  nucleons  are  bound  by 
strong  short-range  forces.  The  forces  can  be  approximated  by  a  deep  rectangular  well. 
The  order  of  the  resulting  levels  is  indicated  iri  the  left  column  of  Fig.  5.  To  arrive  at  the 
empirically  observed  stable  neutron  or  proton  configurations  at  the  magic  numbers  of  2,  8, 
20,  28,  50,  and  82,  a  spin  orbit  coupling  term  was  required;  the  splitting  increases  with  L 
and  the  higher  J-value  lies  deeper.  (See  Ref.  5,  for  example.)  Just  below  the  magic 
numbers  50  and  82,  levels  with  vastly  different  J-values  are  very  close.  Transitions 
between  these  would  lead  to  high  multipole  radiation  and,  as  we  shall  see  below,  these 
have  very  long  lifetimes.  They  are  the  isomers  in  which  we  are  interested.  Years  of 
empirical  observation  have  shown  that,  for  odd  A  nuclei  with  an  odd  number  of  either 
neutrons  or  protons,  there  are  islands  of  isomerism  for  neutron  or  proton  numbers  just 
below  50  and  82.  Isomers  are  also  found  in  other  parts  of  the  periodic  table  for  other 
reasons  (e.g,  shape  isomers). 
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The  electromagnetic  lifetimes  in  nuclei  range  over  many  more  orders  of  magnitude 
than  in  atoms.  A  highly  simplified  nuclear  model  is  used  by  Weisskopf  (Ref.  4)  to  estimate 
the  lifetimes  of  the  various  multipoles  over  a  wide  range  of  energies.  Actual  lifetimes  may 
differ  from  these  estimates  by  one  or  more  orders  of  magnitude.  Figure  6  shows  graphs  of 
these  lifetimes  as  a  function  of  energy  for  both  the  electric  and  magnetic  transitions.  We  are 
interested  mainly  in  the  region  below  100  keV.  Lifetimes  range  from  about  1  fs  to  about 
1020  s. 


Figure  6.  The  half-lives  of  electromagnetic  transitions. 


C.  LINEWIDTHS 


The  width  of  the  transition  is  related  to  the  mean  lifetime  by  the  simple  transform 
relationship: 

T  x  =  h/2:t,  (6) 

where  T  is  the  full  width  at  half  maximum  and  x  is  the  mean  lifetime.  For  a  state  of  high 
excitation,  a  level  may  decay  by  many  paths;  each  path  has  a  partial  width.  Thus: 

r  =  rYi  +  + ...  +  rn  rp  +  Tp  +  ...  (7) 

where  each  subscripted  gamma  is  the  width  for  the  emission  of  that  particular  gamma  ray  or 
particle  and  the  total  width  is  the  sum  of  all  the  partial  widths.  For  a  low  lying  state  in  a 
stable  nucleus,  the  decay  is  generally  entirely  electromagnetic.  It  may  take  place  either  by 
gamma  ray  emission  or  by  internal  conversion  (i.c.).  In  the  case  of  a  low  lying  state,  we 
may  write: 

r*rT+ric  •  (8) 

Internal  conversion  is  a  process  which  may  occur  when  the  excitation  energy  of  the  nucleus 
exceeds  the  binding  energy  of  an  atomic  electron.  In  such  a  case,  the  coulomb  field  may 
mediate  to  transfer  the  excitation  energy  to  the  electron  and  to  expel  that  electron  rather  than 
to  radiate  away  that  energy  by  the  emission  of  a  photon.  This  mechanism  is  the  nuclear 
counterpart  to  the  Auger  process  in  the  atom.  The  conversion  width  is  generally  set 
proportional  to  the  radiation  width  and  is  written  as 

r  =  ry  +  r.  c  =  ry  +  <xr7  =  o  +  a)  rY .  (9) 

The  internal  conversion  coefficient,  a,  is  written  as  a  series  of  coefficients  for  the  various 
shells.  Thus: 

a  =  aK  +  aL  +  aM  +  -  00) 


If  the  excitation  energy  is  between  the  binding  energy  of  the  electron  and  the  mass  of  the 
electron,  there  are  the  following  simplifying  approximations: 


^(E.U-^Z3 


(ID 
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It  is  clear  from  these  approximations  that  the  internal  conversion  coefficient  increases  with 
increasing  atomic  number,  with  increasing  multipolarity,  and  with  decreasing  transition 
energy.  These  effects  show  up  clearly  in  the  more  precise  computations  represented  in  the 
graphs  of  Fig.  7.  The  coefficients  for  the  L  electrons  are  seen  to  have  a  similar  functional 
dependence  on  atomic  number,  multipole,  and  energy.  The  internal  conversion  coefficient 
is  important  because  it  enters  the  formula  for  the  stimulation  cross  section. 


wtf  M.  (Fran  C  M.  L*Mtm  « u,  ~T»b*ot  ha*em.~watt.  Hum  rm.  UMHypniiainJ 


Figure  7.  Internal  conversion  coefficients  for  various  multipole  transitions  in 
nuclei  of  atomic  numbers  30,  60,  and  90. 
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D.  CROSS  SECTIONS 


Consider  transitions  between  the  two  levels  shown  in  Fig.  4.  For  absorption,  the 
cross  section  can  be  obtained  from  the  single  level  Breit-Wigner  formula  and  is  given  by: 


where 


2J2+1  yi  ry 
°abs  2(2Jj  +  1)  It  r 


x  =  (E  -  E  )/(r/2) . 


1 


1+x2 


(12) 


The  stimulation  cross  section  (as  in  the  atomic  case)  is  obtained  from  a  comparison 
of  the  Einstein  B  coefficients.  It  is  given  by: 


2J.  +  1 
°s=  2LM 


abs.  2tc  1+a  i+x2 


=  kX 


1 


(13) 


'2'r*  I+X-  1+X 

A  number  of  factors  affect  the  emission  and  absorption  spectra  of  the  two-level 
system.  These  include  recoil  of  the  absorbing  or  emitting  nucleus  and  Doppler  shifts  of 
frequencies  due  to  the  Maxwellian  distribution  of  velocities  of  the  atoms. 

Recoil  is  important  in  the  nuclear  case  because  the  natural  widths  of  the  transitions 
that  we  shall  consider  will,  in  many  cases,  be  very  narrow.  The  emitting  nucleus  yields  a 
photon  whose  energy  is  the  transition  energy  minus  the  recoil  energy,  Er.  On  the  other 
hand,  the  absorbing  nucleus,  in  its  ground  state,  must  receive  a  photon  whose  energy 
equals  the  transition  energy  plus  the  anticipated  recoil  energy.  Thus,  the  peaks  of  the 
emission  and  absorption  spectra  are  separated  by  twice  the  recoil  energy.  These  conditions 
obtain  for  nuclei  in  the  atoms  of  a  gas  at  0  K.  If  the  gas  is  raised  to  the  temperature  T,  the 
widths  of  the  lines  are  Doppler  broadened.  The  situation  is  depicted  in  the  sketches  of  Fig. 
8.  The  Doppler  width  in  a  gas  is  given  by 
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Figure  8.  Emission  and  absorption  spectra  for  stationary  nuclei.  The  dotted 
curves  represent  the  Doppler-broadened  spectra  for  nuclei  at  the 

temperature  T. 


t 


A  =  E0 


(Ref.  6)  If  the  atoms  are  located  in  a  solid  lattice,  the  expression  is  modified  to: 


(14) 


(15) 


(16) 


where  Cv  is  the  specific  heat  and  ©  is  the  Debye  temperature. 

Generally,  A  »  T ,  and  the  cross  section  for  absorption  (or  stimulation)  becomes: 

°(E)  =  Vf  "a  °°  CXP '  ^rx/2A^  •  <17) 

Recoil  and  broadening  effects  are  compared  for  a  typical  atomic  transition  and  two 
Mossbauer  nuclei.  The  results  are  shown  in  Table  2.  It  is  obvious  from  the  last  column 
that  the  stimulation  cross  section  can  be  diminished  many  orders  of  magnitude  by  Doppler 
broadening. 
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Table  2.  Recoil  and  Broadening  Effects 


^  (cm) 

5>(eV) 

rY  (eV) 

2A  (eV) 

2Er  (eV) 

r/2A 

Na  Atom 

-5 

5.9x10 

2.1 

-8 

4.5  x  10 

(300K) 

6.6  x  10* 

-10 

2x10 

0.007 

191 

Ir 

-10 

9x10 

129k 

-6 

4.7  x  10 

0.14 

0.09 

-5 

3.4  x  10 

57 

Fe 

-9 

8.6  x  10 

14.4k 

-9 

4.6x10 

0.02 

0.004 

-7 

2.3  x  10 

In  order  to  achieve  gain  in  a  pumped  medium,  the  Schawlow-Townes  relationship 
(Ref.  8)  requires  that: 

(N2-g2N1/gi>0]Um>^NV„  ,  (18) 

where  the  N's  and  the  g's  are  the  populations  and  degeneracies  of  the  respective  levels. 
The  ar’s  arc  the  stimulation  and  non-resonant  cross  sections.  The  stimulation  cross  section 
is  given  by: 


sum 


a!  f 


i 


i 


(19) 


2tc  1+a  i+x2  ’ 

where  f  is  the  Mossbauer  fraction,  i.e.,  the  fraction  of  emissions  or  absorptions  for  which 
the  entire  recoil  is  taken  up  by  the  crystalline  lattice.  In  Table  3,  we  compare  the  various 
cross  sections  at  a  range  of  energies.  The  stimulation  cross  sections  exceed  the  non¬ 
resonant  cross  sections;  however,  it  should  be  noted  that  in  a  real  situation,  the  stimulation 
cross  section  would  be  reduced  many  orders  of  magnitude  by  a  combination  of  factors  such 
as  internal  conversion,  Mossbauer  fraction,  Doppler  broadening,  and  crystalline 
inhomogeneous  broadening. 
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Table  3.  Cross  Sections 


E0 

X 

1 

O  , 

photo 

a 

compton 

Pb 

sn 

Pb 

1.2  keV 

10A 

1010b 

■ 

1.5  x  106 

2x10® 

1.5  x  103 

4x103 

12 

i 

00 

o 

■ 

2x10  4 

3x104 

5x102 

2x103 

120 

0.1 

106 

■ 

300 

2x103 

40 

100 

k  =*  [2  jc  (1  +  a)]"1 


E.  STIMULATION  SYSTEMS 

Assume  that,  as  in  Fig.  9,  we  have  a  black  box  containing  an  electromagnetic 
system.  The  box  is  black  only  in  the  visible  so  that  we  cannot  see  into  it.  The  system  is  a 
two-level  system,  capable  of  being  excited  into  the  upper  state  by  absorbing  a  photon  of 
energy  hVa,  and  of  spontaneously  decaying  from  the  upper  state  to  emit  a  photon  of  energy 
hve.  If  the  system  is  in  the  upper  state  and  a  directed  wave  of  energy  hve  passes  over  it,  the 
system  may  be  stimulated  to  emit  a  photon  of  energy  hve  having  the  same  linear  and 
angular  momenta  as  the  stimulating  photon.  What  kind  of  sdmulatable  systems  can  we  find 
to  put  into  the  black  box?  We  shall  describe  two  systems. 

First,  consider  a  two-level  nucleus  in  an  atom  at  rest  If  the  nucleus  is  in  the  upper 
state,  it  can  emit  a  photon  of  energy  (Eq  -  Er),  i.e.,  the  transition  energy  minus  the  recoil 
energy.  The  excited  nucleus  can,  therefore,  also  be  stimulated  to  emit  that  photon  by  a 
directed  photon  of  the  very  same  energy  (Eq  -  Er  ).  However,  a  nucleus  in  the  lower  state 
cannot  absorb  a  photon  of  this  energy;  it  can  only  absorb  a  photon  of  energy  (E0  +  Er).  A 
good  system  is  therefore  a  gas  of  stationary  atoms  with  a  fraction  of  the  nuclei  in  the 
excited  state.  Let  us  consider  only  those  processes  occurring  in  a  direction  parallel  to  the 
main  axis  of  the  acicular  shape.  Any  spontaneously  emitted  photon  can  stimulate  the 
emission  of  another  photon.  A  stationary  nucleus  in  the  lower  level  requires  a  photon  of 
energy  (E0  +  Er)  to  be  excited  to  the  upper  state.  An  atom  whose  nucleus  has  just  been 
stimulated  by  a  stream  of  photons  of  energy  (E0  -  Er)  is  in  the  lower  level  and  is  recoiling 
with  the  energy  Er  towards  the  streaming  photons  with  laboratory  energy  (E0  -  Er).  To 
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absorb  a  photon  and  be  reexcited  into  the  upper  state,  the  nucleus  would  have  to  encounter 
a  photon  of  laboratory  energy  E0.  There  are  none.  There  are  thus  no  nuclear  nonresonant 
processes.  In  this  sense,  the  nuclear  case  differs  from  the  atomic.  Again,  a  viable  system 
is  a  gas  of  atoms  or  molecules  at  rest,  with  each  atom  or  molecule  containing  the  candidate 
nucleus.  Another  system  is  a  gas  whose  constituents  are  all  moving  with  the  same  speed, 
v;  we  need  merely  to  transform  to  a  frame  of  reference  in  which  all  the  atoms  or  molecules 
are  at  rest  and  the  previous  situation  will  again  obtain.  In  other  words,  one  needs  a  system 
whose  constituent  particles  have  a  velocity  distribution  which  is  a  delta  function. 

The  second  example  is  merely  another  way  of  achieving  a  delta  function  velocity 
distribution.  The  atom  is  now  made  part  of  a  strongly  bound  crystalline  lattice.  When  its 
nucleus  emits,  the  recoil  is  taken  up  by  the  entire  crystal,  not  by  the  individual  atom.  The 
recoil  is  therefore  negligible.  In  other  words,  we  require  a  solid  whose  nuclear  emitters 
undergo  a  Mossbauer  effect.  Thus,  both  systems  effectively  have  a  delta  function  velocity 
distribution.  However,  there  is  one  difference.  In  the  gas  case,  the  resonance  energy  for 
stimulation  is  (Eq  -  Er),  and  for  absorption,  (E0  +  Er).  An  excited  state  population  larger 
than  the  ground  state  population  is  not  required  for  gain.  In  the  Mossbauer  case, 
stimulation  and  absorption  occur  at  the  same  energy,  E0;  and,  therefore,  inversion  is 
required. 
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F.  THE  GAS  GRASER 


Two  velocities  immediately  suggest  themselves  -  zero  and  c,  the  speed  of  light. 
The  First,  because  it  is  achieved  by  going  to  0  K;  the  second,  because  it  is  the  upper  limit 
against  which  high-energy  accelerators  squeeze  the  velocity  distribution. 


1 .  The  Velocity  Zero. 


In  1975,  Hansch  and  Schawlow  (Ref.  9)  suggested  that  the  atoms  of  a  gas  could  be 
slowed  down  to  zero  by  irradiating  the  sample  of  gas  from  all  six  directions  with  a  laser 
tuned  to  the  low-frequency  edge  of  a  Doppler-broadened  resonance  line.  Atoms  moving 
towards  the  laser  will  see  its  photons  Doppler-shifted  towards  resonance.  Each  absorption 
will  cause  the  atom  to  lose  an  amount  of  momentum  equal  to  the  momentum  of  the  photon, 
hoc-  The  authors  calculated  that  with  a  2852  A  laser,  a  sample  of  magnesium  gas  could  be 
cooled  from  600  K  to  0.24  K  in  30  jis.  In  1983,  using  only  one  laser  beam  and  a 
counterflowing  beam  of  sodium  atoms,  several  groups  (Refs.  10  and  11)  achieved  a 
stationary  sample  of  gas  with  temperatures  of  the  order  of  50  mK.  Refinements  of  this 
technique  at  the  National  Bureau  of  Standards  (NBS)  and  AT&T  Bell  (Ref.  12)  have 
recently  led  to  samples  of  sodium  gas  at  temperatures  of  about  1  mK. 


Let  us  now  consider  a  gas  sample  of  molecules  (mass  M)  having  a  nuclear 
constituent  with  the  transition  of  interest  characterized  by  the  energy  En  and  the  lifetime  tn- 
The  molecule  also  has  an  atomic  (or  molecular)  transi-ion  which  is  the  cooling  transition;  it 
is  characterized  by  the  transition  energy  and  lifetime  Ec  and  Tc,  respectively.  Let  us 
assume  that  we  are  laser-cooling  from  a  temperature  T  to  0  K.  At  T,  the  mean  velocity  of  a 
molecule  is: 


3  kTl 

l  M 


(20) 


The  decrease  in  velocity  due  to  the  absorption  of  one  photon  of  momentum  AP  is  given  by 

AV  =  AP/M  =  Ec/Mc  .  (21) 

The  number  of  absorptions  to  cool  is  therefore: 

Nc  =  V/AV  =  (3Mc2KT)1/2/Ec  .  (22) 

where  c  is  the  speed  of  light.  One  may  either  pump  the  graser  (or  gamma-ray  laser)  first 
and  then  cool,  or  cool  and  then  pump.  Assume  we  pump  first.  Then,  in  cooling,  assume 
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the  atomic  transition  is  pumped  to  saturation.  Then,  on  the  average,  there  is  one  absorption 
every  Tc.  The  time  to  cool  is  then  given  by  the  number  of  absorptions  times  the  time  for 
one  absorption,  or 

Time  To  Cool  =  Nc  Tc  .  (23) 


To  be  effective  and  not  lose  the  population  inversion,  the  nuclear  lifetime  should  obey  the 
relationship 

%  >  •  (24) 


Given  the  parameters 
T  =  300K, 


M  =  250,  and 


Ec  =  1  eV, 

we  obtain 


For  Tc  =  10'8  s,  we  get: 


tN  >  1.3  x  105  Tc  . 
tN  >  10'3  s  . 


This  situation  could  be  a  severe  constraint  on  the  graser  design.  However,  if  we  cool  first 
and  then  pump,  the  problem  may  be  resolved.  Ultimately,  the  resolution  depends  on  how 
closely  the  temperature  approaches  0  K. 


2.  The  Velocity  c. 

A  high-energy  accelerator  takes  a  particle  from  an  energy  Me2  to  an  energy  yMc2 

where 

o2n-1/2 

Y  =  (1  •  P  )  . 

and 

P  =  V/c  .  (25) 


The  spread  in  energy  is  given  by 

AE  =  (Ay)  Me2  . 


(26) 


The  velocity  spread  is  given  by: 


AP  =  (1/p)  (A y/y3)  =  Ay/y3  (27) 

since  P  ->  1  for  relativistic  particles.  In  the  reference  frame  of  the  beam,  this  expression 
becomes 

Ap'  =  Ay/y  =  AE/E.  (28) 

In  order  to  achieve  a  laser,  the  sample  has  to  be  cooled  to  the  point  where: 

r  ~r 

ad  y 

or 

A3'E0  =  ry .  (29) 

For  a  50  keV  transition  with  Tn  =  10'^  s 

AP’  =  10'12  , 

i.e.. 


AV'  =  0.03  cm/s. 

For  current  accelerators,  the  energy  spread  is  many  orders  of  magnitude  away  from  the 
above  figures.  For  example 


Van  de  Graaffs  [5  MeV]  AP'  =  Ay/y  =  10'5 

Michigan  State  University  Cyclotron  [15  MeV  protons]  10~3 

Bevelac  [1-2  GeV/nucleon]  $  10*2 


(Ref.  13).  Thus,  of  the  two  velocities,  zero  seems  the  more  viable. 

G .  THE  MOSSBAUER  SOLID 

A  group  of  stationary  nuclei  will  display  the  usual  Lorentzian  shape  for  an 

absorption  or  stimulation  resonance,  namely, 

o  =  a  (1  +  x  )  . 

0  v  7 


(30) 


If  the  nuclei  are  in  a  solid,  there  is  a  Maxwellian  distribution  of  velocities,  with  the 
consequent  Doppler  broadening  of  the  resonance.  The  formulae  now  take  the  shape 

a  =  :y-  4  a0  exp  (-^2x2/4)  ,  (31) 

where 

=  r/A  « 1 , 

and 

A  =  E  (2kTeff/Mc2)l/2 , 

and  Teff  is  a  function  of  (©/T)  and  ©  is  the  Debye  temperature.  This  effect,  coupled  with 
recoil,  gives  the  absorption  and  emission  spectra  shown  in  Fig.  10. 


Emission  Absorption 


Figure  10.  Doppler-broadened  emission  and  absorption  spectra. 

The  peaks  of  the  emission  and  absorption  spectra  arc  separated  by  the  energy  interval  2Er. 
At  E0,  there  is  a  sharp  peak  in  each  spectrum  that  is  due  to  the  fraction  of  nuclei  that 
undergo  emission  or  absorption  with  no  individual  recoil  but  have  the  recoil  taken  up  by  the 
entire  lattice.  This  fraction  is  given  by  the  Debye-Waller  factor 

s2  e/T 
1  +1  o'] 

6 


f  =  exp(-2W)  =  exp 


B.J 

'  2k© 


/  _\2  e/T 

•(e)  J  y  <e y  • !) 1  dy 


(32) 


The  fraction  can  approach  unity  for  low-energy  transitions  in  heavy  nuclei  embedded  in 
lattices  with  high  Debye  temperatures.  All  the  current  lasing  schemes  demand  a  Mossbauer 
effect.  The  experimental  Mossbauer  regime  is  approximately  in  the  lifetime  region  from 
10*8  to  about  10-4  s.  That  region  is  being  extended  towards  longer  lifetimes.  The  cross 
section  has  the  form: 


°,tim  =  °0ni  +  x2r 


(33) 
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Sometimes  an  additional  factor  is  applied  to  take  into  account  the  broadening  due  to  the 
inhomogeneities  in  the  crystal. 

H .  PRESENT  APPROACHES  TO  THE  GRASER 

The  three  major  approaches  at  present  are  depicted  in  Fig.  11.  All  begin  with  the 
nuclei  in  a  long-lived  isomeric  state. 

LONG  UPCONVERSION  ELECTRON-NUCLEAR 

LIFETIME  OR  COUPLING 

TWO-STEP 


Isomer 


zGr  zGr  *Gr 

Figure  11.  The  three  major  pumping  schemes. 

In  the  first  approach,  the  isomer  itself  is  used  as  the  upper  state  of  the  lasing  pair. 
Because  of  its  lifetime,  given  that 

PcN  =  6.58  x  10'16  eV-s 

=  0.159  Hz-s,  (34) 

the  radiation  width  is  very  narrow.  In  general,  the  inhomogeneities  of  the  crystal  greatly 
broaden  the  width  and  thereby  reduce  the  effective  cross  section  by  many  orders  of 
magnitude.  One  scheme  is  to  use  nuclear  magnetic  resonance  (NMR)  techniques  to  remove 
the  inhomogeneities.  Other  techniques  are  being  proposed. 

In  the  second  approach,  it  is  assumed  that  there  is  a  level  nearby  which  can  be 
reached  from  the  isomer  by  either  a  laser  or  an  x-ray  source.  This  level  has  a  short  lifetime 
and  cascades  down  to  the  upper  state  of  the  lasing  pair. 

A  two-step  process  is  involved  in  the  third  approach.  First,  one  must  excite  the 
electron  cloud  surrounding  the  nucleus.  Then,  given  compatible  electronic  and  nuclear 
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level  structures,  the  electronic  excitation  transfers  to  the  nucleus.  Again,  either  the  upper 
level  of  the  cascading  pair  has  been  reached,  or  there  is  a  cascade  to  it 

I.  GENERAL  GRASER  PROBLEMS  AND  NEEDS 

If  one  were  to  do  a  systems  study  leading  to  the  actual  design  and  construction  of  a 
graser,  many  studies  would  have  to  be  undertaken  and  many  questions  resolved.  Among 
them  are  the  following: 

1 .  Data  Bases  -  Nuclear  data  bases  specific  to  the  needs  of  the  various  pumping 
schemes  need  to  be  constructed.  Some  are  in  process  or  completed,  (e.g., 
Ref.  14) 

2.  Producing  Macro  Samples  of  Isomers  -  as  experiments  are  proposed, 
samples  of  specific  long-lived  isomers  will  need  to  become  available  in 
sufficient  amounts  to  make  crystals. 

3.  Crystal  Fabrication  -  high-purity,  defect-free  crystals  will  be  needed.  They 
will  have  to  be  very  homogeneous  and  acicular  in  shape.  It  may  be  necessary  to 
implant  the  desired  nuclei  in  a  host  lattice.  Implanting  will  have  to  be 
substitutional.  Crystalline  planes  will  have  to  be  oriented  to  produce  Bragg 
scattering  within  the  crystal  to  enable  the  establishment  of  Borrmann  modes. 

4.  Radiation  Damage  -  the  crystal  will  contain  a  high  density  of  radioactive 
nuclei.  Spontaneous  decay  of  these  nuclei  may  produce  enough  damage  to  the 
lattice  to  decrease  the  Debye-Waller  factor  and  to  produce  inhomogeneous 
broadening.  These  effects  need  to  be  studied. 

5 .  Lattice  Heating  -  in  the  process  of  pumping  the  isomer,  enough  heat  may  be 
generated  in  the  crystal  to  either  destroy  the  Mossbauer  effect  or  to  melt  the 
crystal  itself.  These  effects  are  being  investigated.  (See  Chapter  II.) 

6.  Superradiance  -  many  investigators  assume  that  the  graser  will  be  a 
superradiant  device.  Studies  should  be  conducted  to  determine  whether  the 
conditions  required  for  superradiance  can  be  achieved. 
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J.  SCHEME-SPECIFIC  PROBLEMS 


Each  lasing  scheme  brings  with  it  its  own  characteristic  problems: 

1 .  Long  Lifetime 

•  Mossbauer  Effect  -  Achieving  the  effect  for  long-lived  transitions  has  been 
extremely  difficult.  Experiments  must  be  carried  out,  preferably  using 
candidate  nuclei. 

•  Inhomogeneous  Broadening  -  The  rapid  suppression  of  this  broadening  is 
required  to  trigger  the  energy  release  in  this  scheme.  Demonstration 
experiments  should  be  undertaken  in  the  not  too  distant  future. 

•  Internal  Conversion  -  A  long  lifetime  implies  a  high  multipolarity  which,  in 
turn,  implies  a  large  internal  conversion  coefficient 

2.  Upconversion 

•  Candidate  -  Finding  a  candidate  with  a  suitable  sequence  of  nuclear  levels  may 
present  a  very  difficult  problem.  If  one  cannot  be  found  in  searches  using  the 
current  nuclear  compilations,  theoretical  studies  combined  with  limited 
experimental  searches  will  have  to  be  earned  out 

•  Pumping  -  Laser  pumping  (especially  for  X  >  Ratom  )  may  encounter  extreme 
difficulties  due  to  both  the  impenetrability  of  the  electron  cloud  and  to  lattice 
heating. 

3 .  Electron  •  Nuclear  Excitation 

•  Electron  Excitation  -  The  excitation  of  the  electron  cloud  may  require  the 
deposition  of  sufficient  energy  to  destroy  the  lattice.  Experimental 
demonstrations  are  needed. 

•  Atomic  and  Nuclear  Spectra  -  The  eigenlevel  structure  of  nucleus  and  electrons 
must  be  compatible  to  enable  the  transfer  of  excitation.  Here  again,  a  mixture  of 
theory  and  experiment  may  be  essential. 

•  Mossbauer  Effect  -  The  ability  to  maintain  the  effect  despite  the  deposition  of 
large  amounts  of  energy  during  the  excitation  process  must  be  investigated. 


II.  HEATING  DURING  UPCONVERSION  IN  A 
GAMMA-RAY  LASER 
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A.  INTRODUCTION 

It  has  long  been  hoped  that  very  short  wavelength  lasers  can  be  built  by  using 
nuclear  transitions  to  provide  the  photon  (Ref.  15).  Because  of  the  large  energies  involved, 
it  has  been  proposed  to  store  energy  in  an  isomeric  state,  and  release  it  in  a  single  pulse  by 
upconverting  to  a  lasing  transition  (Ref.  16).  (See  Fig.  1 1  of  this  report.)  It  is  the  purpose 
of  this  note  to  examine  the  heating  of  a  sample  material  during  such  an  upconversion  with 
X-rays.  Heating  is  a  serious  concern  for  a  gamma  ray  laser  because,  unlike  the  situation 
with  X-ray  lasers,  the  Mossbauer  effect  will  be  needed  to  keep  the  photons  and  excited 
states  in  resonance.  The  examination  will  focus  on  the  heating  of  an  57Fe  sample  during 
excitation  of  the  Mossbauer  state  at  14.4  keV.  The  arguments  and  physical  trends  do  not 
extend  to  the  visible  and  near  UV  regime.  These  will  be  examined  in  a  separate  paper. 

To  indicate  how  serious  heating  will  be  during  upconversion  of  nuclear  levels  we 
express  the  law  of  DuLong  and  Petit  for  specific  heat  capacities  of  metals  in  unconventional 
units.  The  standard  value  of  25  J/mole-K  converts  to  2.6x1  O'4  eV/atom-K.  This  means 
that  during  upconversion  with  10  keV  photons,  an  energy  leakage  of  ten  parts  per  million 
going  into  thermal  channels  would  raise  the  temperature  of  the  sample  by  400  K.  To  raise 
iron  from  absolute  zero  to  the  liquid  state  would  require  less  than  65  parts  per  million. 
Resonant  cross  sections  will  be  large,  but  nearly  all  the  absorption  must  be  by  the  nuclear 
transition  in  order  not  to  melt  the  material. 

In  the  second  section  of  this  paper  we  outline  the  calculation,  and  define  the  needed 
parameters  for  the  case  of  ideal  beam— one  which  is  of  zero  width  and  exactly  tuned  to  the 
excitation  energy.  Even  this  favorable  assumption  leads  to  damaging  heating,  except  for 
samples  that  are  very  thin  compared  to  the  radiation  length  of  resonant  photons.  Such  thin 
samples  would  lead  to  inefficient  use  of  the  upconversion  photons  and  increase  pump 
power  requirements. 
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For  non-ideal  triggering  beams  the  situation  is  worse,  because  finite  beam  width 
will  enhance  the  photoelectric  effect  relative  to  the  resonant  excitation.  Expected  heating  as 
a  function  of  beam  width  is  presented  in  Section  C,  along  with  expected  trends  for  variation 
in  upconversion  energy  and  atomic  number.  The  results  indicate  a  serious  problem  which 
would  require  orders  of  magnitude  increases  in  both  beam  intensity  and  resolution. 

It  has  been  suggested  that  these  problems  can  be  dealt  with  by  conductive  cooling, 
(Ref.  17)  or  by  triggering  the  laser  in  a  time  short  compared  to  the  crystal  disintegration 
time  (Ref.  18).  These  ideas  are  discussed  in  Section  D,  along  with  some  preliminary 
results.  Conclusions  are  presented  in  Section  E. 


B.  SAMPLE  CALCULATION— 


In  this  section  we  compute  the  heating  during  upconversion  from  the  ground  state 
to  the  14.4  keV  state  in  57Fe.  We  assume  that  the  incident  radiation  is  exactly  on  resonance 
and  neglect  any  detuning,  and  do  not  include  any  effects  of  heating  due  to  decay  of  the 
excited  state.  For  the  case  calculated,  there  is  no  alternate  decay  channel  except  back  to  the 
ground  state,  and  a  relatively  small  internal  conversion  coefficient.  This  approach  results  in 
a  very  favorable  case  from  the  point  of  view  of  heating. 


The  essential  ingredients  in  this  calculation  are  the  resonance  cross  section  for 
upconversion  of  the  desired  state,  the  competing  photoelectric  cross  section,  along  with  the 
energy  loss  cross  sections  for  the  resulting  electrons,  and  finally  the  geometry  of  the 
sample.  Knowledge  of  the  basic  bulk  properties  of  iron  is  also  needed.  First  we  assume 
that  the  nuclear  excitation  is  exactly  on  resonance.  The  cross  section  for  excitation-de¬ 
excitation  is: 


rer 


(rT) 


2  ’ 


(35) 


where  X  is  the  wavelength,  and  the  Fs  are  the  elastic  (Te),  inelastic  (Tj)  and  total  widths 
(Pr).  For  the  case  under  consideration  here,  simple  upconversion  to  a  state  which  only 
decays  back  to  the  ground  state,  the  cross  section  simplifies  to: 

2  re 


(36) 


The  ratio  of  widths  is  often  written  1/(1 +a),  where  a  is  the  internal  conversion  coefficient. 
This  is  the  ratio  of  internal  conversion  to  photon  emission  decays  of  the  state.  For  the  case 
under  consideration,  the  resonant  cross  section  is  2.38  Mb  (Ref.  19).  The  competing 
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photoelectric  cross  section  is  only  about  7000  b,  (Ref.  20)  or  200  times  smaller.  The 
corresponding  radiation  lengths  (Lx)  defined  as  the  reciprocal  of  the  linear  extinction 
coefficient,  are  given  in  Table  4. 


Table  4.  Radiation  lengths  associated  with  resonance  excitation,  photoelectric 
effect  and  high  energy  (14  keV)  electron  loss  in  57Fe. 


57  Fe 

l/L(A‘l) 

L  (pm) 

°R 

1.11E-03 

0.0899 

aPE 

5.94E-06 

16.8 

(dE/dx)e 

7.90E-05 

1.27 

To  compute  the  heating  associated  with  a  given  level  of  upconversion  it  is  necessary 
to  know  the  probability  of  photoelectron  emission  as  a  function  of  the  probability  of 
l  resonance  excitation.  Because  excitation  of  the  resonant  state  removes  potential  targets 

from  future  photons,  the  relative  effectiveness  of  the  photoelectric  process  increases  as  the 
population  is  inverted.  The  correct  expression  for  the  total  probability  of  photoelectron 
(Ppe)  emission  as  function  of  the  final  fraction  of  nuclei  in  the  resonant  state  (Pr)  is 


L« 

pPE^ta 


1-P. 


This  in  an  important  correction  for  Pr  =  0.5  or  greater. 


The  next  physical  parameter  that  is  needed  is  the  average  energy  deposited  in  the 
material  per  photoelectron.  This  energy  comes  from  two  sources.  First,  the  photoelectron 
itself  deposits  energy  on  its  way  out.  Second,  the  residual  excited  atom  decays  either  by 
electromagnetic  or  Auger  processes,  and  these  electrons  and  photons  deposit  energy.  In 
order  to  ensure  that  the  heating  was  not  being  overestimated,  we  assumed  that  the  photo¬ 
effect  produced  a  high-energy  electron  and  left  relatively  little  energy  behind  in  the  excited 
atom.  Since  energy  loss  per  unit  length  increases  rapidly  as  the  electron  slows  down, 
assuming  a  high  energy  photoelectron  underestimates  the  calculated  energy  deposition, 
provided  the  sample  is  small  enough  that  the  photoelectrons  escape. 
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This  assumption  of  an  energetic  electron  also  leads  to  a  minimal  contribution  to 
heating  from  the  atomic  excitation,  since  relatively  little  energy  is  deposited.  To 
characterize  this  numerically  we  have  taken  the  effective  energy  absorption  cross  section  to 
be  1Mb,  which  is  the  cross  section  of  100  eV  photons  in  carbon  (Ref.  21).  We  expect  this 
assumption  to  underestimate  the  true  absorption.  To  turn  this  number  into  a  rough  energy 
deposition,  we  compute  Eatom  =  Ppe  E*[l  -  exp  (-  a/2L)]  where  E*  is  the  excitation 
energy  of  the  atom  (the  incident  x-ray  energy  minus  the  ejected  electron's  energy)  and  a  is 
the  radius  of  the  wire.  The  energy  deposition  per  atom  estimated  in  this  fashion  is  given  in 
Table  5. 

The  energetic  photoelectron  will  have  a  range  in  the  material  which  varies  as  the 
energy  squared.  In  other  words,  it  loses  energy  progressively  more  rapidly  as  it  slows 
down.  This  is  observed  empirically  and  the  trend  persists  down  to  energies  on  the  order  of 
100  eV.  We  obtain  a  lower  limit  on  the  energy  deposition  by  linearizing  the  energy  loss. 
For  photoelectrons  uniformly  produced  in  a  long  cylinder,  with  isotropic  velocities,  the 
angle  averaged  distance  to  travel  to  the  edge  of  the  cylinder,  s,  increases  from  the  radius  for 
particles  at  the  edge  to  jc/2  times  the  radius  for  particles  at  the  center.  We  can  compute  a 
lower  limit  on  the  energy  deposited  by  the  photoelectron  by  using  the  minimum  value  of  s, 
along  with  a  radiation  length  associated  with  the  initial  energy  of  the  particle  to  write 

kpE  =  PPE  df  -  Smin  * 

Since  photoelectrons  are  preferentially  produced  perpendicular  to  the  direction  of  motion  of 
the  incident  photon,  the  assumption  of  isotropy  further  decreases  the  estimated  energy 
deposition.  We  hope  to  remove  these  approximations  in  order  to  get  a  more  reliable 
estimate.  At  the  moment  it  is  worth  noting  that  our  present  "estimate"  is  in  fact  a  very 
conservative  lower  bound.  The  energy  depositions  from  the  energetic  electrons  are  also 
given  in  Table  5  for  a  variety  of  iron  wire  thicknesses,  along  with  the  total  values  from  both 
contributions.  Also  shown  is  the  final  temperature  of  the  sample,  assuming  an  initial 
temperature  of  100  K,  and  a  heat  capacity  from  the  Law  of  DuLong  and  Petit  This  serves 
to  make  the  numbers  roughly  applicable  for  any  metal. 

Now  let  us  consider  the  consequences  of  the  results  in  Table  5.  For  a  variety  of 
radii  of  iron  wire  (first  line)  we  have  computed  the  effective  absorption  probability  (second 
line).  All  calculations  assume  an  incident  photon  number  equal  to  the  number  of  atoms  in 
the  sample.  The  fraction  of  the  sample  upconverted  increases  with  thickness,  but 
somewhat  slower  than  linearly.  The  fraction  of  photons  which  convert  to  photoelectrons 
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increases  linearly.  The  photoelectric  energy  deposition  is  approximately  quadratic  in  the 
thickness.  The  final  temperature  increases  dramatically  with  sample  thickness  as  shown  in 
Fig.  12.  A  cold  (100  K)  sample  of  thickness  0.01  pm  would  not  heat  up  to  room 
temperature;  however,  doubling  the  thickness  to  0.02  pm  results  in  a  final  temperature 


Table  5.  Heating  as  a  function  of  whisker  thickness  for  5?Fe.  In  aii  calculations 
the  fluence  is  scaled  so  that  the  photon  number  equals  the  number  of  atoms  in 

the  sample. 


a  (Hm) 

0.01 

0.0201 

0.032 

0.036 

0.0375 

0.04 

0.05 

0.0625 

Pr 

0.11 

0.20 

0.30 

0.33 

0.34 

0.36 

0.43 

0.50 

PPE 

0.0006 

0.0012 

0.0019 

0.0021 

0.0022 

0.0024 

0.0030 

0.0037 

DEPOSITION 
FROM  PHOTO¬ 
ELECTRON  (eV) 

0.03 

0.13 

0.34 

0.43 

0.47 

0.53 

0.83 

1.29 

DEPOSITION 
FROM  ATOM 
(eV) 

0.01 

0.06 

0.15 

0.19 

0.20 

0.23 

0.36 

0.57 

TOTAL 

0.05 

0.19 

0.49 

0.62 

0.67 

0.76 

1.19 

1.86 

TF 

283 

840 

1810 

1836 

1945 

2137 

3030 

4425 

E-OUT  (mJ) 

0.02 

0.16 

0.60 

0.84 

0.984 

1.12 

2.08 

3.82 

Figure  12.  Temperature  increase  with  whisker  thickness  for  perfectly 

on*resonance  beam 
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of  840  K,  which  is  almost  twice  the  Debye  temperature.  By  0.03  (im  the  sample  has  begun 
to  melt. 

There  are  two  difficulties  with  avoiding  heating  by  going  to  very  small  samples. 
First,  one  would  have  difficulties  in  preparation.  Second,  the  energy  stored  becomes  very 
small.  It  is  certainly  true  that  using  nuclear  energy  levels  allows  very  efficient  storage  of 
energy,  but  if  only  a  small  fraction  of  a  small  sample  can  be  triggered,  then  there  might  be 
no  practical  applications.  Table  5  also  gives  the  released  energy,  assuming  the  given 
conversion  efficiency.  We  see  that  very  small  energies  are  involved,  small  even  compared 
to  a  single  pulse  from  a  conventional  pulsed  laser. 

C.  UPCONVERSION  WITH  FINITE  WIDTH  BEAMS 

The  effect  of  finite  beam  width  on  the  upconversion  process  is  to  make  the  effective 
resonance  cross  section  smaller,  while  leaving  the  photoelectric  process  essentially 
unchanged.  Synchrotron  radiation  may  be  used  to  provide  well  resolved  X-rays  in  the  1  to 
10  keV  energy  range,  with  beam  widths  on  the  order  of  0.1  eV,  and  intensities  of  5  x  108 
photons  per  burst.  These  pulses  would  typically  last  around  20  ns,  with  a  pulse  repetition 
frequency  of  1000  Hz. 

In  order  to  effectively  upconvert  a  nuclear  state,  one  needs  a  lifetime  long  compared 
to  the  length  of  the  burst,  and  an  intrinsic  width  large  compared  to  the  beam  spread.  Since 
a  10  ns  lifetime  corresponds  to  a  6.6  x  10-8  eV  intrinsic  width,  and  these  quantities  have  a 
reciprocal  relationship  it  is  impossible  to  simultaneously  meet  both  requirements  with 
current  beams.  Recently,  a  very  high  resolution  (0.005  eV)  measurement  was  reported, 
(Ref.  22)  but  this  resolution  improvement  came  at  the  cost  of  greatly  reduced  intensity,  as 
is  usually  the  case.  The  most  noteworthy  features  of  these  beams  are  that  they  are  broad 
compared  to  the  width  of  the  state,  they  are  not  intense  enough,  and  finally  there  are 
relatively  few  photons  within  one  natural  linewidth. 

In  the  case  of  upconversion  with  a  broad  beam,  it  is  necessary  to  account  for  the 
decreased  effective  resonance  cross  section.  In  general,  it  is  also  necessary  to  account  for 
"bleaching"  of  the  beam,  i.e.,  preferential  removal  of  photons  from  the  resonance  peak,  but 
for  very  thin  samples  as  considered  here,  this  effect  can  be  neglected.  For  large  beam 
spread  AE  compared  to  the  width  of  the  state,  the  average  cross  section  varies  like 
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The  1/AE  dependence  is  due  to  the  factor  normalizing  the  beam  to  fixed  intensity.  The 
constant  of  proportionality  depends  on  the  functional  form  used  to  represent  the  beam 
shape.  Differences  in  this  constant  are  small.  We  used  a  uniform  energy  distribution  of 
total  width  AE,  which  gives  for  the  effective  cross  section 
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The  effect  of  finite  beam  width  is  shown  in  Fig.  13,  where  the  energy  deposition  in 
57Fe  is  plotted  as  a  function  of  the  ratio  of  beam  spread  to  natural  width.  The  target  wire  is 
assumed  to  have  a  0.01  pm  diameter.  The  upconversion  was  set  at  1  percent  of  the 
sample.  This  is  the  minimum  consistent  with  gain  during  stimulated  emission.  The  arrow 
indicates  the  ratio  appropriate  for  a  0.1  eV  width  and  a  state  with  a  10  ns  lifetime.  It  is  clear 
that  there  is  too  much  energy  deposition  by  orders  of  magnitude. 

It  has  long  been  recognized  that  the  Borrmann  effect  would  be  needed  to  reduce  the 
heating  through  atomic  scattering.  An  optimistic  assessment  might  anticipate  a  factor  of 
100  suppression  of  the  photoelectric  cross  section.  The  calculations  were  repeated 
assuming  such  a  suppression,  and  the  results  are  presented  in  Fig.  14.  Even  under  these 
circumstances,  with  current  resolution,  nearly  100  times  the  energy  necessary  to  melt  the 
sample  would  be  deposited  in  the  process  of  upconverting  1  percent  of  the  nuclei. 

The  energy  of  the  calculation  was  chosen  from  the  observed  transition.  We  can 
ask,  what  would  the  heating  be  for  a  lower  energy  X-ray  transition?  Naively,  one  might 
expect  decreased  energy  deposition,  but  in  fact  the  opposite  results.  With  softer  X-rays, 
the  photoelectric  cross  section  increases  dramatically,  generally  scaling  as  1/E4.  The 
resonance  cross  section,  meanwhile,  will  not  change  much  at  the  peak  because  of  internal 
conversion.  The  factor  a  from  the  ratio  iyrr  in  equation  (36)  increases  like  A2,  hence  the 
peak  cross  section  is  essentially  constant.  To  make  matters  worse,  the  width  of  the  state 
decreases  as  the  energy  of  the  transition  decreases.  This  means  that  for  fixed  resolving 
power,  the  ratio  of  AE/T  will  increase  with  lower  energy.  (This  change  in  width  is  not 
included  in  the  calculations  presented  here,  since  results  are  shown  as  a  function  of  AE/T.) 
In  addition  to  lower  energy  x-rays  producing  more  photoelectrons,  the  resulting  low  energy 
electrons  lose  their  energy  faster  on  the  way  out  of  the  material. 

Expected  heating  of  an  iron-like  material  by  1  keV  upconverting  photons  is  shown 
in  Fig.  15  for  the  same  target  whisker  as  above.  A  beam  width  equal  to  a  few  times  the 
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Figure  13.  Energy  deposition  in  eV  per  atom  during  1%  upconversion  of  a 
0.01  iim  57Fe  sample.  Current  beam  resolution  is  inadequate  to  prevent 
severe  overheating,  even  with  beams  of  insufficient  intensity. 
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Figure  14.  Energy  deposition  in  eV  per  atom.  Same  conditions  as 
Fig.  13,  but  with  a  factor  of  100  reduction  in  the  photoelectric  cross  section. 

Overheating  is  still  a  problem. 


3. 


Figure  15.  Same  as  Fig.  14,  but  for  an  iron-like  target  with  a  1  keV  transition 
with  the  same  nuclear  matrix  elements  as  the  14.4  keV  transition.  The  heating 
has  become  much  more  severe  because  the  photoelectric  effect,  and  the  rate  at 
which  the  photoeiectrons  deposit  energy  are  both  dramatically  increased. 


natural  width  of  the  state  already  contributes  enough  energy  to  heat  the  sample  above  the 
Debye  temperature. 

This  situation  will  continue  to  get  worse  to  energies  below  100  eV.  It  may  begin  to 
improve  at  that  point  because  dE/dx  of  tne  photoeiectrons  begins  to  decrease. 
Nevertheless,  the  heating  is  probably  intolerable  for  any  photons  above  the  minimum 
ionization  threshold. 

D.  LOOPHOLES 

In  this  section  we  discuss  two  proposed  loopholes,  conductive  cooling,  and  time 
considerations  which  might  allow  triggering  of  a  gamma-ray  laser,  despite  the  heating. 
Conductive  cooling  of  the  sample  requires  that  the  sample  be  embedded  in  or  against  a 
material  which  has  a  high  thermal  conductivity,  high  heat  capacity,  and  low  photoelectric 
cross  section.  We  considered  carbon  as  a  candidate  heat  sink.  (Other  materials  will  be 
considered  in  the  future.)  We  took  the  Debye  temperature  as  a  target  temperature,  and 
computed  the  energy  needed  to  raise  graphite  from  around  100  K  to  440  K,  the  Debye 
temperature  of  iron.  Then,  for  given  AE/r  we  computed  the  energy  deposited  per  iron 
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atom  (assuming  the  Borrmann  effect,  as  in  Fig.  14),  and  converted  that  to  the  number  of 
carbon  atoms  needed  per  iron  atom,  and  a  thickness  for  the  carbon  backing.  We  then 
computed  the  direct  heating  of  the  carbon  backing  due  to  the  photoelectric  effect  in  the 
backing.  The  carbon  was  also  assumed  to  be  in  a  crystal  with  a  factor  of  100  reduction  in 
the  photoelectric  cross  section  due  to  the  Borrmann  effect.  No  account  was  taken  of 
photoelectrons  from  one  material  passing  into  and  heating  the  other. 

Despite  these  optimistic  assumptions,  the  results  shown  in  Table  6  indicate  limited, 
if  any,  help  from  the  heat  sink.  For  0. 1  eV  resolution  beam,  we  find  that  the  needed 
carbon  heat  sink  is  sufficiently  large  that  despite  the  lower  photoelectric  effect,  it  would  be 
heated  above  the  temperature  of  the  iron,  and  aggravate  the  problem.  As  the  postulated 
beam  is  improved,  the  carbon  begins  to  lower  the  temperature  of  the  iron.  Nevertheless, 
there  appears  to  be  no  appreciable  region  where  this  approach  would  be  helpful.  It  should 
be  stressed  that  this  again  is  a  conservative  calculation,  and  that  substantial  heating  of  the 
carbon  from  photoelectrons  produced  in  iron,  and  vice  versa  should  be  expected.  On  the 
other  hand,  berylium  or  boron,  because  of  their  even  smaller  photoelectric  cross  sections 
might  be  better  in  this  application  as  heat  sinks.  Future  calculations  will  address  this 
question. 


Table  6.  Heatsink  temperature  (K)  due  to  photoelectric  effect  in  the  heatsink 
itself.  Size  of  heatsink  set  by  excess  heating  (above  Tqebye)  in  0.01  pm  Fe 

sample,  initially  at  100  K. 


ae/t 

Ed 

T  (K) 

a  (pm) 

Ed(C) 

T  (C) 

106 

23.0 

0.32 

mm 

Hotter 

105 

2.3 

5346 

0.1 

ESI 

Cooler 

104 

0.23 

985 

0.032 

320* 

103 

0.023 

189 

0.01 

mu 

♦Insufficient  for  cooling  of  entire  sample  below  Tdebye- 

Another  approach  to  evading  the  difficulties  presented  by  heating  is  to  attempt  to  get 
the  coherent  radiation  out  of  the  lasing  medium  before  the  heat  leads  to  disintegration. 

Generally,  the  lifetimes  under  consideration  are  on  the  order  of  10  ns.  If  they  are 
much  shorter,  then  pump  power  problems  are  believed  to  become  insurmountable.  The 
question  then  is,  will  the  crystal  hold  together  well  enough  for  both  the  Mossbauer  and 


Borrmann  effects  to  apply  for  a  time  on  the  order  of  10  ns.  We  believe  that  this  is  at 
present  an  unanswered  question.  From  the  results  of  Ref.  23  it  is  clear  that  the  time 
required  for  break  up  of  a  diffraction  pattern  is  a  very  strong  function  of  the  amount  of 
energy  deposited.  The  data  are  plotted,  along  with  a  simple  power  law  fall-off,  in  Fig.  16. 
(The  best  fit  was  t  =  6723  x  ^  ps,  where  x  is  the  energy  absorbed  divided  by  the  energy 
required  to  melt  the  sample.)  This  time  decreases  very  rapidly  with  deposition  energy.  If 
barely  enough  energy  is  deposited  for  melting  of  the  sample,  the  crystal  integrity  times  may 
be  on  the  order  of  a  few  ns  in  aluminum.  If  three  times  the  melting  energy  is  deposited,  the 
crystal  integrity  time  extrapolates  to  a  few  ps.  Since  diffraction  depends  on  fixed  position, 
whereas  the  Mossbauer  effect  relies  on  low  velocity,  the  Mossbauer  effect  would  probably 
be  lost  on  a  time  scale  short  compared  to  that  for  the  diffraction  pattern.  Characteristic 
times  for  electron  transit,  even  at  low  energy,  are  quite  short  compared  to  the  likely 
lifetimes.  (A  100  eV  electron  travels  at  6  (X/ps.)  The  combination  of  crystal  disintegration 
times  and  electron  transit  times  both  short  compared  to  the  relevant  nuclear  lifetimes  makes 
lasing  before  disintegration  appear  improbable.  Further  analysis  will  include  electron- 
photon  collision  times  and  atomic  excited  state  photon  and  auger  decay  times. 


ENERGY  ABSORBED/ENERGY  TO  MELT 


Figure  16.  Crystal  disintegration  times  from  Ref.  23  as  a  function  of  deposited 
energy.  The  energy  deposition  is  in  units  of  the  energy  required  to  melt  the 

sample. 
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E.  CONCLUSIONS 

X-ray  upconversion  of  a  nuclear  isomer  to  a  suitable  lasing  state  is  going  to  be  an 
extremely  difficult  undertaking.  Intrinsically,  leakage  of  a  very  small  fraction  of  energy 
into  thermal  channels,  primarily  via  the  photoelectric  effect,  will  destructively  heat  the 
sample.  Current  X-ray  sources  are  orders  of  magnitude  too  broad  in  energy,  as  well  as  too 
weak  to  accomplish  the  task.  It  appears  that  the  time  scales  of  the  destructive  heating  are 
short  compared  to  nuclear  decay  times.  A  preliminary  examination  indicates  that 
conductive  cooling  of  the  sample  is  not  a  workable  option.  Photoelectric  heating  of  the 
heatsink  itself  prevents  its  utility. 

From  the  above  discussion  it  appears  that  two  courses  of  action  are  called  for. 
First,  within  the  context  of  the  upconversion  approach  to  building  a  gamma-ray  laser,  much 
more  attention  should  be  paid  to  possible  triggering  beams.  Furthermore,  other  approaches 
not  relying  on  upconversion  should  be  receiving  greater  attention. 
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III.  AN  INVESTIGATION  OF  THE  DESTRUCTION  BY 
INHOMOGENEOUS  BROADENING  OF  RESONANCE  IN 
ISOMERIC  CRYSTALS  AND  ITS  RESTORATION  BY 

SPECIAL  EFFECTS 

A.  INTRODUCTION 

The  first  concepts  for  developing  a  gamma-ray  laser  were  based  on  stimulated 
emission  from  isomeric  crystals  (Ref.  24).  It  was  natural  to  examine  the  possibility  of 
stimulated  emission  in  long-lived  nuclear  transitions  after  (1)  the  development  of  the  maser 
and  laser  showed  that  inversion  was  possible  and  stimulated  emission  could  be  used  to 
generate  and  control  coherent  microwave  and  optical  radiation  (Ref.  6)  and  (2)  the 
discovery  of  the  Mossbauer  effect  showed  that  nuclear  resonance  fluorescence  could  be 
easily  observed  under  well-defined  conditions  (Ref.  25).  The  possibility  of  using  nuclear 
isomeric  crystals  as  the  lasing  media,  the  isomeric  level  being  the  upper  lasing  level,  was 
desirable  because  it  would  allow  the  use  of  lower  pumping  powers  for  the  inversion. 

However,  soon  after  the  discovery  of  the  Mossbauer  effect,  experimental 
exploration  of  this  phenomenon  of  nuclear  resonance  revealed  that  observing  the  effect  with 
nuclei  of  longer  lifetimes  (narrower  natural  linewidths)  become  increasingly  more  difficult 
(Ref.  13).  The  reason  for  this  is  that  nuclei  at  different  sites  in  the  crystal  experienced 
slightly  different  electric  and  magnetic  fields  because  of  crystal  irregularities,  impurities, 
and  even  different  nuclear  spin  orientations  of  nearest  neighbors.  These  effects, 
collectively  referred  to  as  inhomogeneous  broadening  of  the  Mossbauer  line,  will  be 
discussed  in  more  detail  later.  For  the  present,  it  is  sufficient  to  point  out  that,  in  order  to 
use  isomeric  levels  in  a  lasing  transition,  the  inhomogeneous  broadening  effect  on  the 
resonance  has  to  be  reduced.  Two  concepts  have  been  proposed  to  do  just  this. 

One  concept  for  producing  a  gamma-ray  laser  is  based  on  the  ability  to  remove  the 
effect  of  the  coupling  which  causes  inhomogeneously  broadened  nuclear  resonance  lines  in 
the  isomeric  crystal.  The  technique  proposed  for  eliminating  the  effects  of  nuclear  dipole- 


dipole  interaction,  a  major  contributor  to  the  inhomogeneous  broadening,  is  a  well-known 
RF  pulsing  technique  used  in  high-resolution  nuclear  magnetic  resonance  (NMR) 
spectroscopy  (Ref.  26)  and  theoretically  shown  to  be  applicable  to  Mossbauer  spectroscopy 
(Refs.  27  and  28).  Such  techniques  have  been  used  to  reduce  the  inhomogeneous 
broadening  in  NMR  work  by  up  to  four  orders  of  magnitude  (Ref.  29). 

In  another  concept,  the  individual  narrow  lines  are  homogeneously  broadened. 
Overlap  of  lines  and  resonance  is  achieved  (Ref.  30)  because  as  a  result  of  this  operation  all 
the  nuclei  in  the  sample  have  the  same  broad  lines,  which  are  greater  than  the 
inhomogeneous  broadening. 

Level  Mixing  Spectroscopy  (LEMS)  is  based  on  the  principle  of  homogeneous  line 
broadening  (Ref.  31).  Recently,  the  application  of  this  technique  to  a  perturbed  angular 
correlation  (PAC)  type  of  experiment  led  to  a  measurement  of  a  "resonance"  in  109Ag.* 
The  effect  of  this  level  mixing  on  a  Mossbauer  type  of  experiment  has  been  discussed  in 
recent  publications  (Refs.  30,  and  32).  Experiments  are  being  planned  to  verify  these 
concepts  (personal  communication,  P.  Boolchand). 

Another  way  of  obtaining  homogeneous  broadening  of  resonance  lines  which  leads 
to  greater  overlap  between  source  emissions  and  the  absorber  cross  section  is  through 
relaxation  phenomena  or  time-dependent  hyperfine  interactions  (Ref.  33  and  34).  These 
homogeneous  broadening  effects  are  well  known  in  the  Mossbauer  community.  Their  use 
in  recovering  the  resonance  destroyed  by  inhomogeneous  broadening  in  the  source  and 
absorber  experiment  has  to  be  investigated. 

In  this  paper  we  examine  the  utility  of  line  shifting  and  broadening  techniques  to 
increase  the  resonance  effect  in  inhomogeneously  broadened  systems. 

B.  RESONANT  LINESHAPES  AND  MOSSBAUER  EXPERIMENTS 

The  interaction  between  a  source  and  an  absorber  in  a  Mossbauer  experiment  can  be 
described  in  terms  of  an  emission  line  and  an  absorption  cross  section.  If  a  detector  is 
located  as  shown  in  Fig.  17,  then  for  a  thin  absorber  which  does  not  appreciably  change 
the  beam  line  profile,  and  assuming  only  resonant  absorption, 
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Strictly  speaking  this  is  not  a  resonance  in  the  Mossbauer  sense. 
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Figure  17. 


Mdssbauer  transmission  experimental  geometry. 


the  detected  loss  in  the  total  photon  count,  1^,  because  of  the  presence  of  the  absorber  is 
given  by 

h.  +•» 

ID=  J  dt  J  I  (E,r,S)  n  a(E,r)  dE  (37) 

t,  - 

where  n  is  the  number  of  resonant  nuclei  in  the  beam  per  cm2,  and  the  measurement  is 
performed  during  the  time  interval  tj  -  tj.  I  (EX.S)  is  the  source  radiation  lineshape  and 
o(E,D  is  the  resonant  cross  section  of  the  nuclei  in  the  absorber. 


The  source  lineshape  is  given  by 
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i  (E,r,s)  =  X  pi  (ei> 
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[E-(E0-ei)-S]2  +  (r/2)2 


(38) 


where  Iq  is  the  flux  of  photons  from  the  source  in  cm'2  s'1  in  the  direction  of  the  detector, 
•  E0  is  the  resonant  energy  of  an  unperturbed  nucleus,  £j  is  the  perturbation  of  the  energy  of 

nucleus  i  due  to  the  environment,  S  is  the  Doppler  shift  due  to  the  relative  motion  of  the 
source  and  absorber  nuclei  and  is  given  by  S  =  v/c  Eq,  with  velocity  v  where  c  is  the  speed 
of  light,  T  is  the  natural  linewidth*  and  fs  is  the  recoil-free  fraction  in  the  source.  It  is 


* 


T  is  the  total  level  width  which  is  a  sum  of  the  level  widths  of  all  the  contributing  processes  to  the 
level  decay. 
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assumed  in  equation  (37)  that  the  integration  time  is  short  enough  that  IQ  can  be  taken  as 
constant. 

The  distribution  of  nuclear  environments  is  assumed  to  be  normal  with  width  A  so 
that  the  probability  of  finding  a  nucleus  with  energy  e*  displaced  from  the  resonance  energy 
E0  in  a  width  Aq  is  given  by 


1 


P  (e.)  =  • 

11  /ft  A 


-eV 


Ae. 


Similarly,  for  the  effective  cross  section  of  the  absorber  we  have 

a0rj-fa 

o(E,T)  =  X  Pi  <ei>  - ~'i — ~f~T 

i  [E-(E0-ej)]2  +  (y)2 


(39) 


(40) 


Where  c0  is  the  cross  section  on  resonance  and  is  given  by 

1 


a 

0  2k  21+1  1  +  a  ’ 


g 

Ie  and  Ig  are  the  nuclear  spin  quantum  numbers  of  the  excited  and  ground  states, 
respectively,  a  is  the  internal  conversion  coefficient  of  the  transition,  X  is  the  nominal 
wavelength  of  the  gamma-ray,  Ty  is  the  photon  absorption  partial  width.  Figure  18  is  a 
schematic  representation  of  the  resonance  effect  when  Ei  =  o.  The  source  and  absorber 
nuclear  resonance  energies  are  given  by  Es  and  Ea,  respectively,  and  AE  =  Ea  -  Es.  The 
integral  over  energy  given  in  equation  (37)  is  represented  by  the  shaded  region  in  Fig.  18. 


Figure  18. 


Resonance  condition  for  source  and  absorber 
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The  usual  way  of  performing  Mossbauer  experiments  is  depicted  in  Fig.  17.  One 
moves  the  source  at  predetermined  velocities  and  counts  the  number  of  transmitted 
y-photons.  In  this  way  a  Mossbauer  transmission  spectrum  is  generated  to  obtain 
information  about  the  nuclear  electromagnetic  interactions. 

The  Mossbauer  effect  is  defined  in  (Ref.  35)  as 

R(S  =  qq)-R(S) 

(S)  R  (S  =  oo) 

with  the  observed  counting  rate  for  a  "thick"  absorber  given  by 

oo 

R  (S)  =  J I  (E,r,S)  exp  [-J!t  (E)  T|  dE,  (41) 

o 

where 

pt  (E)  =  n fa  (E,H  +  |ie  • 

In  equation  (41)  |ie  is  the  linear  extinction  coefficient  which  is  a  function  of  the 
photoelectric  effect  and  Compton  scattering  cross  sections,  n  is  the  number  of  resonant 
nuclei  per  cm2,  T  is  the  thickness  of  the  absorber  in  cm,  (3  =  fna0T  is  usually  referred  to  as 
the  equivalent  absorption  thickness  or  the  Mossbauer  thickness  of  the  absorber  (Ref  35). 

There  are  other  ways  of  moving  the  resonances  apart  in  a  controlled  fashion  than 
that  shown  in  Fig.  17.  One  can  introduce  (1)  external  DC  magnetic  fields  to  shift  the 
resonances;  (2)  vary  the  temperature  of  the  samples  and  thus  affect  changes  in  fa,  fs,  and  to 
a  much  smaller  extent,  E0  (second-order  Doppler  shift);  or  (3)  use  RF  fields  to  produce 
zero-order  shifts  in  the  energy  of  the  nuclear  levels  and  mix  the  nuclear  states  (produce 
dressed  states  of  nuclei),  as  suggested  by  Collins  et  al.  (Ref.  36). 

C.  INHOMOGENEOUS  BROADENING 

If  the  absorber  and  source  were  composed  of  nuclei  subjected  to  the  same  internal 
and  external  fields  at  the  same  temperature,  the  result  of  the  experiment  would  be 
essentially  as  described  in  equation  (37)  with  Ej  =  0.  The  signal  strength  would  be  a 
function  of  the  number  of  nuclei  in  the  path  of  the  beam,  pa  T  A,  where  pa  is  the  nuclear 
density  in  the  absorber,  T  is  the  thickness  of  the  absorber,  and  A  the  area  of  the  absorber 
illuminated  by  the  beam.  This  assumes  no  multiple  scattering  and  a  thin  sample,  so  that 
there  is  no  modification  of  the  beam  as  it  traverses  the  absorber  (effects  already  well  studied 
by  other  investigators)  (Refs.  35,  37).  The  actual  experimental  conditions  are  generally 
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quite  different.  Nuclei  imbedded  in  a  solid  host  interact  with  the  neighboring  atoms  and 
change  this  ideal  situation  drastically.  The  effect  of  impurities,  crystal  dislocations,  and 
crystal  boundaries  is  to  introduce  slightly  different  fields  at  nuclei  located  in  different 
positions  in  the  lattice  and  thereby  shift  the  nuclear  resonant  energies.  The  effects  are  most 
pronounced  for  isomeric  (long-lived)  levels. 

Even  if  the  crystal  were  perfect  and  we  considered  only  those  nuclei  well  within  the 
regular  structure,  so  that  boundary  effects  were  negligible,  there  would  still  be  a  slightly 
different  resonance  energy  at  different  lattice  sites.  The  nuclear  spin-spin  or  magnetic 
dipole-dipole  interaction  between  the  emitting  or  absorbing  nuclei  and  their  neighbors 
would  be  different  due  to  random  orientations  of  the  magnetic  moments  of  different  nuclei. 
This  is  an  inherently  nuclear  effect  and  if  it  were  possible  to  eliminate  all  the  solid  state  and 
geometrical  factors,  this  interaction  would  still  remain. 

The  spin-spin  interaction  energy  is  given  by 


r 


where  Pi,  M-2  are  the  magnetic  moments  of  the  respective  nuclei  and  r  is  the  distance 
separating  their  centers.  For  107 Ag,  this  energy  has  been  estimated  to  be  about  10" 12  eV 
(Ref.  38).  Such  a  weak  interaction  is  not  important  in  experiments  with  good  Mossbauer 
nuclei,  such  as  57Fe,  which  has  a  natural  linewidth  of  5  x  10'9  eV.  It  has  been  argued  that 
this  interaction  essentially  destroys  resonance  between  long-lived  or  isomeric  nuclei  such  as 
107  Ag  or  109Ag  because  the  effective  lineshape  for  the  absorber  and  source  is  now  smeared 
out  over  energies  much  larger  than  the  natural  linewidth,  which  is  on  the  order  of  10'17  eV. 

Inhomogeneous  broadening  of  resonance  levels  is  depicted  in  Fig.  19,  where  the 
resonance  condition  for  a  good  Mossbauer  isotope  57Fe  is  compared  to  the  resonance 
condition  for  a  poor  Mossbauer  isotope, 107 Ag.*  In  view  (a)  of  Fig.  19,  the  resonances  of 
the  emitting  and  absorbing  nuclei  are  separated  by  an  energy,  Ae,  which  arises  because  of 
the  slight  difference  in  the  local  environment  of  the  two  nuclei.  The  natural  linewidth,  F  of 
57Fe,  is  large  enough  that  the  overlap  region  provides  a  good  resonance  effect.  On  the 
other  hand,  as  shown  in  view  (b)  of  the  figure,  for  the  long-lived  isotope  with  the  very 
narrow  natural  line,  the  environmental  difference  between  the  nuclear  sites,  which  induces 
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For  the  moment,  the  difference  in  the  recoilless  fraction,  f,  between  the  two  isotopes  is  ignored. 


the  energy  shift  Ae,  destroys  the  resonance  condition  between  emitter  and  absorber.  To 
properly  describe  this  situation,  £j  *  0  has  to  be  introduced  into  the  equations. 


Figure  19.  Resonance  conditions  for  short-  and  long-lived  isotopes.  For 
comparison  we  show  in  (a)  57  f©  a  good  Mfissbauer  isotope  but  difficult  to  invert 
because  of  the  short  lifetime  and  in  (b)  107Ag,  relatively  more  easy  to  invert 
because  of  the  long  lifetime  but  not  a  good  Mdssbauer  isotope  because  of  the 

narrow  lines. 

This  parameter  can  be  conveniently  divided  into  three  parts  so  that 


+  en.  +  Ee. 
1  1 


where  £q  is  the  energy  shift  in  the  resonance  of  nucleus  i  due  to  crystal  imperfections, 
impurities,  boundary  effects,  eni  is  the  energy  shift  in  the  resonance  of  nucleus  i  due  to 
•  nuclear  spin  effects  from  neighboring  nuclei  and,  Eei  is  due  to  external  effects  such  as 

gravitational  differences  at  different  nuclei  sites. 

The  average  effect  on  Mossbauer  experiments  of  this  inhomogeneous  line 
broadening  due  to  the  energy  shifts  £;  can  be  approximated  by  an  effective  increase  in  the 
®  linewidth  of  the  source  or  absorber  by  a  factor  (1  +  a)  where  a  is  called  the 

"inhomogeneous  broadening  parameter"  (Ref.  13).  The  effect  on  the  threshold  condition 
for  lasing  is  drastic,  because  the  factor  (1  +  a )  multiplies  the  total  linewidth  of  the  nuclei  in 
the  lasing  medium,  thus  effectively  reducing  the  resonance  cross  section. 
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In  the  line-narrowing  concept  for  developing  a  7- ray  laser  (Ref.  39),  the  parameter  a 
has  to  be  substantially  reduced.  It  is  assumed  that  the  line-broadening  terms  due  to  crystal 
imperfections  and  impurities  can  be  reduced  by  proper  geometry,  size,  and  crystal 
prepartion  and  that  the  nuclear  spin  interaction  term  can  be  reduced  by  an  appropriate 
sequence  of  RF  pulses,  as  in  high-energy  resolution  NMR.  Alternatively,  a  technique  for 
homogeneously  broadening  the  individual  resonance  lines  such  as  Level  Mixing 
Spectroscopy  (LEMS)  can  be  applied  to  reduce  a.  The  two  mechanisms  for  achieving 
resonance  in  the  sample,  or  more  precisely  getting  more  nuclei  on  resonance  are  depicted  in 
Fig.  20.  In  case  one  the  narrow  lines  are  shifted  toward  E,  the  unperturbed  energy,  and  in 
case  two  the  lines  are  broadened  so  that  the  effective  linewidth  is  greater  than  the  energy 
producing  the  inhomogeneous  broadening.  The  total  magnitude  of  the  required  effect  is 
dependent  on  other  isomer  characteristics,  as  discussed  in  Ref.  (39). 


Figure  20.  Two  different  approaches  for  obtaining  resonance  with  long- 
lived  isotopes  (destruction  of  the  effect  of  inhomogeneous  broadening)  in  (a) 
an  RF  pulsing  technique  causes  a  line  shift  and  in  (b)  through  cross  fields  level 
mixing,  or  relaxation  homogeneous  broadening  of  lines  causes  resonance 

overlap. 

The  inhomogeneous  broadening  term  a  will  be  reduced  naturally  when  either  the 
ground  or  excited  nuclear  states  have  zero  spin,  because  then  either  pi  or  (12  will  be  zero. 
Also,  if  the  nuclei  in  the  crystal  are  all  initially  in  the  excited  state  the  effect  of  this  term  will 
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be  greatly  diminished  in  the  early  stages  of  laser  pulse  evolution  because  of  the  relatively 
small  number  of  nuclei  in  the  ground  state. 


D .  THE  MOSSBAUER  RESONANCE  IN  THE  PRESENCE  OF 
INHOMOGENEOUS  BROADENING 


Assuming  the  source  lineshape,  absorber  cross  section,  and  distribution  of  nuclear 
environments  as  given  by  equations  (38,  39,  and  40,  respectively)  the  intensity  of  the  beam 
I’  absorbed  by  a  thin  absorber*  (E0  =  0  without  loss  of  generality)  is  given  by 

T  H  N. 

I'(T,A,S)  =  —2-2-  n,  £ 


X 


J 


dE 


AkA^A^ 


(I72tc) 


e  e  Ae.  Ae. 


(E-Ej-S)2  4-  (T/2)2  (E-e.)2  +  (r/2)2 


(42) 


where  Ns,  Na  are  the  total  number  of  sites  with  energy  shifts  ej  and  £j  in  the  source  and 
absorber,  respectively,  and  where  na  is  the  number  of  nuclei  per  cm2  in  the  absorber. 
Converting  the  sum  to  an  integral  one  has  (see  Appendix  A) 


I’(r,A,S)  =  &  IA0T/A)  na  Re  (er)  erfc(Z) 


(43) 


where 


A2  s  (A2  +  A2)/2  ,  and  Z  s  (T  +  iS)//2  A  . 

The  maximum  normalized  resonant  absorption  (S  =  0)  is  given  by 

I"(r, A)  =  I’(r,A,S  =  0)/  ( j  I0nacr0) 

whercx=7f7 ' 


(44) 


*  A  thin  MOssbauer  absorber  is  one  for  which  the  lineshape  of  the  beam  passing  through  the  absorber  is 
not  sufficiently  modified  to  show  saturation  effects.  This  is  discussed  in  Appendix  B,  where  it  is 
shown  that  for  the  MOssbauer  thickness  b  <  10  the  ratios  of  the  intensities  of  the  lines  in  the  spectrum 
are  close  to  their  theoretical  values  and  saturation  effects  are  small.  Under  such  conditions,  it  is 
appropriate  to  use  equation  (37)  for  the  absorption  calculations  instead  of  the  general  expression  given 
in  equation  (41). 


47 


The  definite  integrals  presented  in  Appendix  C  are  helpful  in  examining  some  special  cases 
of  equations  (42)  and  (43).  We  now  examine  some  special  limiting  cases  assuming  Ai  = 

^2- 

Case  (a),  rare  packing.  For  all  i  and  j  in  the  source  and  absorber,  T<<  lej  -  ej  I. 

If,  in  addition  to  this  assumption,  the  number  of  nuclei  in  the  sample  is  small  such 
that,  on  the  average,  the  near  resonance  distance  is  nT,  where  n  is  the  number  of  linewidths 
between  resonances  as  shown  in  Fig.  20,  equation  (42)  can  be  reduced  to: 


^o^ona 

f  r  ^ 

ly 

2A27C 

rY 

2n2 

r  ‘ 

2  2 
-£* /A 


Case  b,  dense  packing  with  T  »  lei  -  Ej  I,  for  all  i. 

The  linewidth  is  much  greater  than  the  separation  between  the  resonances.  This 
implies  that  T  =  A.  Then, 


Case  c,  dense  packing  with  T,  S  «  A 


r  (r,A)  = 


m 

4 


Iq^0^3 
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Case  d,  dense  packing  with  T  or  S  »  A 

L  t  I"., 


i'(r,A) « 


y 

XT2  +  S2) 


The  normalized  maximum  resonance  absorption,  I"(r,A)  calculated  from  equation  (44)  is 
plotted  in  Fig.  21  as  a  function  of  17 A  for  the  case  T/Ty  =  1,  and  in  Fig.  22  as  a  function  of 
A/T  for  various  values  of  T/Fy.  For  inhomogeneous  broadening  on  the  order  of  the  natural 
linewidth,  T/A  =  1,  good  absorption  is  obtained,  indicating  the  presence  of  a  strong 
Mossbauer  resonance  between  the  source  and  absorber.  For  A  »  T,  the  absorption  drops 
linearly. 


9-21-87-11 


Figure  21.  A  plot  of  the  maximum  normalized  resonance  absorption  I" 
(calculated  from  equation  (44),  assuming  Ty/r  =  1)  as  a  function  of  r/A  . 
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NORMALIZED  RESONANCE  ABSORPTION  INTENSITY,  l,,(r,A) 


Figure  22.  Plots  of  the  maximum  normalized  resonance  absorption  curves,  • 

nr, A),  calculated  from  equation  (44)  and  plotted  as  a  function  of  A/r.  The 
numbers  labeling  the  curves  give  the  different  values  of  ryr.  The  dashed 
horizontal  arrow  shows  the  effect  on  the  resonance  expected  from  RF  pulsing 
techniques  (Refs.  26  and  29)  and  the  solid  arrow  the  reduction  in 
inhomogeneous  broadening  expected  from  good  crystal  preparation  techniques 

from  ordinary  isomeric  samples  with  a  broadening  parameter  of  a  =  106.  • 
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E.  DECREASING  THE  EFFECT  OF  INHOMOGENEOUS  BROADENING 
IN  MOSSBAUER  EXPERIMENTS 

Inhomogeneous  broadening  results  from  nuclear  environments  at  different  lattice  sites 
producing  different  shifts  in  the  resonance  energies  £i  of  the  respective  nuclei.  This  effect 
reduces  the  overlap  in  the  source  and  absorber  lineshapes.  In  order  to  increase  the  overlap, 
one  can  either  shift  the  individual  resonances  back  to  a  common  energy  E0,  as  shown  in 
Fig.  20(a)  or  broaden  the  lines  homogeneously  without  shifting,  as  shown  in  Fig.  20(b). 
It  may  be  possible  to  perform  a  modification  of  the  lineshape  which  results  in  both  effects 
as  shown  in  20(c).  Two  approaches  to  overcome  the  inhomogeneous  broadening  effect  on 
Mossbauer  experiments  have  been  presented  in  the  literature.  Generally,  the  first  step  in 
both  techniques  is  the  preparation  of  a  crystal  of  the  isomer,  with 
minimum  deformations.  However,  even  in  an  ideal  crystal  the  dipole-dipole  interaction 
between  non-zero  spin  nuclei  will  persist  and  produce  different  energy  shifts  at  different 
nuclear  positions. 

In  the  first  approach,  one  attempts  to  reduce  the  effect  of  the  dipole-dipole 
interaction  by  the  technique  of  pulsed  RF  fields.  This  technique  was  originally  developed 
for  high-resolution  NMR  spectrometry  and  is  incorporated  into  commercial  instrumentation 
(Ref.  26).  It  has  been  claimed  to  reduce  inhomogeneous  linewidth  in  nuclear  ground  state 
spin  transitions  by  as  much  as  104  (Ref.  29)  and  has  been  examined  theoretically  as  a 
procedure  for  reducing  the  inhomogeneous  broadening  in  Mossbauer  effect  experiments 
(Ref.  27).  It  reduces  the  inhomogeneous  linewidth  by  breaking  the  effect  of  the  dipole- 
dipole  interaction  through  rapid  reorientation  of  the  spins,  thus  effectively  shifting  lines 
toward  the  same  resonance  energy  Ae— »Ae'— >0  for  all  i)  as  shown  in  Fig.  20(a).  With 
reference  to  Fig.  22,  the  resulting  decrease  in  A  would  move  the  operating  point  from  point 
(1)  to  point  (2)  on  the  curve,  reflecting  a  decrease  in  A/T  of  104.  To  approach  the  desired 
operating  point  (3)  a  further  reduction  in  A  would  have  to  be  achieved  by  proper 
preparation  for  the  crystal  to  minimize  the  effect  of  impurities  and  crystal  imperfections. 

The  second  technique  is  depicted  in  Fig.  20(b).  Instead  of  shifting  the  resonance 
energy  at  each  nucleus  by  -£i  to  coincide  with  the  base  nuclear  resonant  energy  [see 
equation  (42)]  the  second  technique  homogeneously  broadens  each  line  to  obtain  greater 
overlap  in  the  lineshapes  as  shown  in  Fig.  20(b).  The  homogeneous  line  broadening  can 
be  achieved  by  various  means.  In  particular  by  working  in  the  region  of  level  crossings, 
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mixing  of  sublevels  can  provide  enhanced  transitions  between  sublevels  leading  to 
homogeneous  broadening  of  lines. 


The  reduction  of  the  overlap  of  nuclear  resonance  lines  in  an  inhomogeneously 
broadened  crystal  is  depicted  in  Fig.  23(a).  Only  some  of  the  lines  are  shown  for  clarity 
presentation.  Partial  recovery  of  the  overlap  due  to  line  shifts,  homogenous  broadening, 
and  a  combination  of  the  two  effects  is  depicted  in  Fig.  23(b)  through  23(g). 

Level  mixing  effects  have  been  used  in  Level  Mixing  Spectroscopy  (LEMS)  to 
obtain  resonances  in  perturbed  angular  correlation  experiments  in  109 Ag  (Ref.  31).  Such 
resonances  have  not  been  seen  before  in  the  Ag  transition  and  are  not  seen  outside  of  the 
level  crossing  region,  indicating  the  importance  of  the  shortened  coherence  time  in  these 
experiments.  In  Mossbauer  experiments,  on  the  other  hand,  the  value  of  the  resonance 
cross  section  determines  the  size  of  the  effects.  If  this  broadening  is  done  at  the  expense  of 
the  maximum  cross  section  on  resonance,  then  as  can  be  seen  from  equations  (40)  and  (43) 
an  increase  in  the  total  width  T  will  not  be  effective  in  producing  an  increase  in  resonance 
as  T  -4  A  because  the  factor  Ty/T  in  equation  (44)  will  still  dominate  even  though  17 A  the 
argument  in  the  exponential  and  error  function  factors  increases  to  values  greater  than  one  . 
The  homogeneous  broadening  produces  overlap  of  lines  but  this  is  not  sufficient  to  get 
good  absorption  because  of  the  factor  Ty/T.  The  mechanism  (LEMS)  has  to  induce  level 
shifts  in  addition  to  homogeneous  broadening  in  order  to  be  effective. 

Another  way  of  obtaining  homogeneous  broadening  is  through  relaxation  effects 
which  are  clearly  exhibited  in  Mossbauer  effect  spectra.  Figure  24  shows  Mossbauer 
transmission  spectra  of  Ferrichrome  A,  a  paramagnetic  compound,  at  several  temperatures. 
The  spectra  range  from  a  well-defined,  six-line  pattern  with  narrow  lines  approaching  the 
natural  iinewidth  at  low  temperatures  to  a  complex  broadened  pattern  with  linewidths 
determined  by  the  electron  spin  relaxation  time  and  the  magnitude  of  the  interaction  to  a 
single  motionally  narrowed  peak,  with  the  hyperfine  pattern  destroyed  at  high 
temperatures. 
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RELATIVE  INTENSITY 


DISTANCE  FROM  UNPERTURBED  RESONANCE 
IN  UNITS  OF  THE  NATURAL  LINE  WIDTH 


Figure  23.  Reduction  of  the  overlap  of  nuclear  lineshapes  in  an 
9  inhomogeneous  crystal  (a)  and  partial  recovery  of  the  overlap  due  to  line  shifts 

(b,  c)  homogeneous  broadening  (d,  e)  and  a  combination  of  the  two  (f,  g).  Note 
the  decrease  in  the  relative  intensities  as  the  width  Increases  (d,e,f,g). 
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Figure  24.  Nuclear  lineshape  modification  due  to  electronic  relaxation 
without  change  in  the  nuclear  lifetime:  experimental  results  in  Ferrichrome  A 
are  shown,  together  with  the  theoretical  calculation.  Source  of  data  is  Ref.  34. 
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Figure  25  shows  a  hypothetical  source  absorber  overlap  of  lines  with  narrow  lines 
in  Fig.  25(a)  and  broadened  lines  in  Fig.  25(b),  and  25(c).  For  the  case  depicted  in  Fig. 
25(a)  the  Mossbauer  effect,  which  is  dependent  on  the  overlap  integral,  is  proportional  to 
rYr/AE2,  whereas  in  cases  (b)  and  (c)  it  is  proportioned  to  (TyAE);  an  increase  by  a  factor 
of  AE/r. 
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Figure  25.  Nuclear  lineshape  modification  due  to  electronic  relaxation 
without  change  in  the  nuclear  lifetime  for  two  different  relaxation  roles.  The 
solid  and  dashed  curves  represent  the  source  and  absorber  lineshapes, 
respectively.  The  shaded  region  Indicates  the  overlap. 
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F.  CONCLUSIONS 


In  a  source-absorber  experiment  where  both  the  source  and  absorber  are 
inhomogeneously  broadened  and  the  distribution  of  fields  at  the  individual  nuclear  sites  is 
normal,  the  Mossbauer  effect  is  inversely  proportional  to  the  width  of  the  distribution.  The 
resonance  effect  can  be  restored,  in  general,  by  removing  the  energy  shifts  at  the  individual 
sites  but  not  simply  by  homogeneously  broadening  the  lines,  even  though  significant 
overlap  is  achieved.  Under  special  conditions,  homogeneous  broadening  of  lines  will 
provide  an  increase  in  the  Mossbauer  effect.  Such  special  conditions  can  occur  if  the 
broadening  is  accompanied  by  a  significant  lineshift  toward  the  mean  of  the  line  positions 
or  if  the  broadening  produces  an  interaction  between  nuclei,  leading  to  cooperative  effects 
and  a  speed  up  of  the  radiative  decay,  as  in  Dicke  superradiance. 
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APPENDIX  A 


EVALUATION  OF  THE  THIN  ABSORBER  INTENSITY 


In  Section  DI-D,  the  intensity  of  a  beam,  I’,  transmitted  by  a  thin  absorber  was 
considered;  the  results  were  used  primarily  for  the  special  case  of  S=0,  that  is,  coincident 
distributions  of  lines  in  the  source  and  absorber.  In  this  appendix,  the  details  of  the 
derivation  of  the  beam  intensity  are  given  for  all  S  and  for  absorber  and  source  line 
distributions  of  different  widths.  In  integral  form,  the  intensity  for  unit  density  (na  =  1)  is 
given  by 
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where  A2  and  A2  are  the  source  and  absorber  widths,  respectively.  The  integral  over  E  is 
easily  performed  by  contour  integration. 
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The  remaining  integrations  are  simplified  by  a  change  of  variables  to: 
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In  these  variables,  the  x  and  y  integrals  can  be  done  independently. 
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where  2A2  =  Ai2  +  A22.  The  integral  over  x  is  trivial,  leaving  the  following  integral  for 
the  intensity: 
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This  final  integral  is  performed  by  rewriting  the  denominator  via  a  Laplace  transform: 
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where  Re  denotes  the  real  part  of  the  integral.  The  integration  over  y  is  now  trivial  and  the 
final  answer  is  given  in  terms  of  the  integral  over  s  which  is  easily  recognized  as  the 
complementary  error  function.  Thus,  the  intensity  of  the  beam  is  given  by: 

r  =  m  Re  erfc(z) 
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This  exact  combination  of  exponential  and  error  function  is  tabulated  in  standard  references 
(for  example,  "Handbook  of  Mathematical  Functions",  edited  by  Abramowitz  and  Stegun; 
note  that  in  this  reference  z  is  replaced  by  -iz,  effectively  interchanging  the  real  and 
imaginary  parts).  | 


The  expression  given  above  can  be  confirmed  in  several  limits.  For  A  =>  0 
(T  and/or  S  =>  °°),  the  gaussians  in  ei  and  E2  in  the  original  form  of  the  integral  can  be 
replaced  with  delta  functions  giving  the  simple  expression 
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which  is  also  obtained  from  the  above  expression  if  it  is  noted  that  for  large  z 

ez2  erfc  (z)  =  — - —  r  (A. 9) 

Jk  z 


I 


A-4 


and  taking  the  real  part  On  the  other  hand,  for  A  =>  (r  and  S  =>  0 ),  the  exponential  in 
the  final  integral  over  y  can  be  neglected,  yielding 

r  =  (a.  io) 

A 

in  agreement  with  the  above  expression  for  z  =  0. 
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A  DISCUSSION  OF  THE  EFFECTS  ON  LINESHAPE  OF  THE 
MOSSBAUER  THICKNESS  PARAMETER 

The  Mossbauer  thickness  parameter  is  defined  as  (5  =  nOofT  where  n  is  the  number 
of  nuclei  per  cm3  in  the  beam,  o0  is  the  maximum  resonant  cross  sections,  f  is  the 
recoilless  fraction,  and  T  is  the  thickness  of  the  sample.  The  Mossbauer  spectrum  for  large 
(3  is  affected  by  the  thickness  in  a  complex  way  due  to  saturation  effects  but  generally  for  p 
<  10  such  effects  are  unimportant.  Figure  B-l(a)  shows  a  Mossbauer  scattering  spectrum 
obtained  from  an  iron  bar  with  a  natural  content  of  57Fe  (=2  percent).  The  ratio  of  line 
intensities  approaches  the  3:2:1  ratio  expected  theoretically  for  small  p.  In  Fig.  B-l(b),  for 
comparison,  the  results  for  a  sample  enriched  in  57Fe  is  shown.  The  saturation  effects  are 
clearly  evident  when  comparison  is  made  with  the  expected  line  intensities  shown  in  Fig. 
B-2.  In  both  Figs.  B-l(a)  and  B-l(b)  the  solid  curves  represent  calculated  spectra  for  the 
respective  cases  (Ref.  B-l). 


FIGURE  B-1.  (a)  Experimental  and  theoretical  scattering  results  for  a 
1/8-in. -thick  iron  bar  at  room  temperature  having  the  natural  content  of  Fe 
(2%).  The  dots  are  the  data  and  the  solid  curve  is  the  calculated  result.  The 
only  free  parameter  is  the  percent  effect  at  one  peak,  (b)  Same  kind  of  r®su*ts_ 
as  shown  above  but  the  scatterer  is  a  90%  enriched  iron  powder  ( f  =  0.8  and  p  = 

175)  at  room  temperature  (Ref.  B-1). 
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DEFINITE  INTEGRALS 
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